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Context. Knowing how the mass loss of carbon-rich AGB stars depends on stellar parameters is crucial for stellar evolution modelling, 
as well as for the understanding of when and how circumstellar structures emerge around these stars, e.g., dust shells and so-called 
detached shells of expelled gas. 

Aims. The purpose of this paper is to explore the stellar parameter space using a numerical radiation hydrodynamic (RHD) model 
of carbon-star atmospheres, including a detailed description of dust formation and frequency-dependent radiative transfer, in order to 
determine how the mass loss of carbon stars changes with stellar parameters. 

Methods. We have computed a grid of 900 numeric dynamic model atmospheres (DMAs) using a well-tested computer code. This 
grid of models covers most of the expected combinations of stellar parameters, which are made up of the stellar temperature, the 
stellar luminosity, the stellar mass, the abundance of condensible carbon, and the velocity amplitude of the pulsation. 
Results. The resultant mass-loss rates and wind speeds are clearly affected by the choice of stellar temperature, mass, luminosity and 
the abundance of available carbon. In certain parts of the parameter space there is also an inevitable mass-loss threshold, below which 
a dust-driven wind is not possible. Contrary to some previous studies, we find a strong dependence on the abundance of free carbon, 
which turns out to be a critical parameter Furthermore, we have found that the dust grains that form in the atmosphere may grow too 
large for the commonly used small particle approximation of the dust opacity to be strictly valid. This may have some bearing on the 
wind properties, although further study of this problem is needed before quantitative conclusions can be drawn. 
Conclusions. The wind properties show relatively simple dependences on stellar parameters above the mass-loss threshold, while the 
threshold itself is of a more complicated nature. Hence, we chose not to derive any simplistic mass-loss formula, but rather provide a 
mass-loss prescription in the form of an easy-to-use FORTRAN routine (available at http://coolstars.astro.uu.se ). Since this mass-loss 
routine is based on data coming from an essentially self-consistent model of mass loss, it may therefore serve as a better mass-loss 
prescription for stellar evolution calculations than empirical formulae. Furthermore, we conclude that there are still some issues that 
need to be investigated, such as the role of grain-sizes. 
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1. Introduction 

How the mass loss of carbon-rich AGB stars (carbon stars) de- 
pends on stellar parameters is not very well known. It is cru- 
cial, however, to have that kind of information in many con- 
texts, such as for stellar evolution modelling. A variety of 



empirical (e.g., \^ssiliadis & Wood 19931 : iGro enewegen 1993': 
Groenewegen 1995b as well as theoretical (e.g.. Blocker 1 995.: 
Arndtetal. 19 97: Wac hter et al. 20021: IWachter etal 2003) for- 



mulae with different parameterisations are available in the liter- 
ature and the agreement between them is rather poor. These sim- 
ple mass-loss prescriptions are nevertheless used in many mod- 
els of stellar evolution, sometimes without enough consideration 
of their applicability in specific cases. 

Determining the mass-loss rates observationally is compli- 
cated and riddled with many possible sources of uncertainty. 
Also, the fundamental stellar parameters, such as effective tem- 
perature and mass, can be hard to determine observationally with 
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high precision, which makes it very difficult to empirically deter- 
mine the mass-loss rate as a function of stellar parameters with 
sufficient accuracy. Empirical studies as such can only provide 
a limited amount of information, since one cannot easily dis- 
entangle the effects of individual stellar parameters. The reason 
for this is of course that stellar evolution and observed stellar pa- 
rameters are connected, and we can only observe each star at one 
evolutionary stage. Observations can tell us about the mass loss 
for a given set of stellar parameters, but the exact dependences 
on these parameters can unfortunately not be determined easily 
due to selection effects. Existing empirical relations suffer from 
selection effects, with a very steep dependence of mass loss rates 
on stellar parameters. 

The fundamental idea behind a so-called dust driven wind 
is that stellar photons, incident on dust grains, will lead to a 
transfer of momentum from the radiation to the atmospheric gas, 
which is dragged along by the dust grains. Models of this process 
have existed for quite some time and have become increasingly 
more complex. Wood (1979) and Bowen (1988) made dynamic 
atmosphere models introducing a parameterised opacity to de- 
scribe the effects of dust formation in the circumstellar envelope. 
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Bowen & Willson (1991) further elaborated on the implications 
of their models for the mass loss during AGB evolution, i.e., the 
development of a so-called "superwind" phase. Blocker (1995) 
then derived a modified version of Reimers' law based on the 
models of Bowen (1988) to prescribe mass loss in his calcula- 
tions of stellar evolution on the AGB. This mass-loss prescrip- 
tion is still widely used, although it is entirely based on the early 
work by Bowen (1988). 

Fleischer et al. (1992) presented self-consistent dynamic 
models of carbon rich AGB stars that include a time-dependent 
description of dust formation. Their work revealed phenomena 
like a discrete spatial structure of the circumstellar dust shell and 
multiperiodicity. A picture has emerged with dust-driven winds 
as a complex phenomenon, where the onset of mass loss requires 
particular circumstances. For about a decade, time-dependent 
dynamical models for stellar winds (e.g., Hofner & Dorfi 1997; 
IWinters et al. 2000) were based on rather crude descriptions of 
radiative transfer and simplified micro-physics, which led to un- 
realistic atmospheric density-temperature structures. 

Hofner et al. (2003) presented dynamic model atmo- 
spheres which couple time-dependent dynamics and frequency- 
dependent radiative transfer, including a non-grey descrip- 
tion of the dust component. Comparing the observable prop- 
erties calculated from these state-of-the-art dynamical mod- 
els (spectra and their variations over time) with observa- 
tions has shown drastic improvements compared to the syn- 
thetic spectr a based on previous (grey ) generations of mod- 



els (see, e.g., Gautschv-Loidl et al. 2004t iNowotnv et al. 2005at 



iNowotnv et al. 2005bl) . While further improvement may become 
possible in the future, a reasonable level of realism has now been 
reached in the modelling to allow application to stellar evolution. 
The model used in the present study is identical to that used in 
Mattsson et al. (2007a) and contains only minor modifications of 
the Hofner et al. (2003) model. It is probably the most advanced 
theoretical tool available at present for studying how dust driven 
winds of carbon stars depend on stellar parameters. 

Detailed numerical modelling provides valuable insights into 
the mass loss problem, since it can be used to constrain the ac- 
tual physics involved in wind formation. In this paper, we ex- 
plore the stellar parameter space using a numerical radiation hy- 
drodynamic (RHD) model of carbon-star atmospheres, including 
time-dependent dust formation and frequency-dependent radia- 
tive transfer, in order to determine how the mass-loss of carbon 
stars depends on stellar parameters. The results we give here re- 
flect primarily the intense mass loss during the so-called thermal 
pulse phase, i.e., the very late stages of carbon-rich AGB evolu- 
tion. 



2. Theory and methods 

2.1. RHD atmosphere models 

The model includes frequency-dependent radiative transfer and 
non-equilibrium dust formation, i.e., we solve the coupled sys- 
tem of frequency-dependent radiation hydrodynamics (RHD) 
and time-dependent dust formation employing an implicit nu- 
merical method and an adaptive grid (cf. Hofner et al. 2003). 

In our model, the stellar atmosphere and circumstellar enve- 
lope are described in terms of conservation laws for the gas, the 
dust and the radiation field, expressed by the following set of 
coupled, nonlinear partial differential equations (PDEs): 

- The three equations describing conservation of mass, mo- 
mentum and energy for the gas. 



- The 0th and 1 st moment equation of radiative transfer. 

- Four moment equations of dust formation (cf. Section |23] l. 

- The Poisson equation (self gravity). 

To our system of nonlinear partial differential equations de- 
scribing the physics of a C-star atmosphere we add a so called 
"grid equation" which determines the locations of the grid points 
according to accuracy considerations (Dorfi & Drury 1987) and 
an equation keeping track of the amount of condensible carbon, 
leaving us with a total of 12 partial differential equations (PDEs). 
This system of PDEs is then solved implicitly using a Newton- 
Raphson scheme. All equations are discretised in a volume- 
integrated conservation form on a staggered mesh. The spatial 
discretisation of the advection term is a monotonic second-order 
advection scheme (van Leer 1977). The same order of numerical 
precision is used for all PDEs. Details of the numerical method 
are discussed by Dorfi & Feuchtinger (1995) and in several pre- 
vious papers about dust-driven wind models (cf. Hofner et al. 
1995 and references therein). 

2.2. General procedure 

The RHD computations are started from hydrostatic dust-free 
initial models providing the atmospheric structure at f = 0. When 
the dust equations are switched on, dust condensation starts and 
the resulting radiative acceleration creates an outward motion of 
the dust and the gas. In the first computational phase, the expan- 
sion of the atmospheric layers is followed by the grid to about 
20 - 30/?* (usually around to ~ 10'^ cm). At this distance, the 
location of the outer boundary is fixed, allowing outflow. The 
outflow model then evolves typically for more than 100 peri- 
ods. To avoid a significant depletion of mass inside the com- 
putational domain, the model calculation is stopped after about 
10^ time steps. This is important since we cannot (due to the 
computational method) allow for material to flow over the inner 
boundary. Due to this fact, a large enough mass loss will even- 
tually lead to a depletion of mass in the considered part of the 
atmosphere and circumstellar envelope. 

2.3. Inner boundary: Stellar pulsations 

The dynamic atmosphere model does not include a physical 
model for the pulsation mechanis m. The p ulsation is modelled 
by a "piston boundary condition" jBowen 1988,) located at Rm ~ 
0.8 - Q.9Ri, in order to be in the optically thick part of the atmo- 
sphere but outside the zone where the pulsations are excited. It 
should be noted that all results obtained in this paper are, to some 
extent, subject to this inner boundary condition, and as shown 
by Mattsson et al. (2008), the effect on the mass-loss rate is not 
negligible in the non-saturated wind regime where the degree of 
dust condensation is affected by the input of kinetic energy by 
pulsations. 

We use a harmonic piston boundary given by the expressions 



R,n(t) ^Rin(0) + — sin (ojt)., 
a> 

and 



Min(f) — AUp COS (cot) , 



2n 



(1) 



(2) 



where !P is the period and Aup is the velocity amplitude 
of the pulsations. To restrict the number of free parame- 
ters, we also employ an empirical period-luminosity relation 
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dFeast et al 1989h and keep the period V tied to the luminosity 
for all models. The piston amplitude Au^, or the kinetic-energy 
injection by pulsation (cf. the ^-parameter in Mattsson et al. 
2008) 



1 



(3) 



is essentiall y a free parameter in t he model, but previou s expe- 
rience (e.g. jHofner & Dorfi 1997t iMattsson et al. 2007ah shows 
that only a limited range can be considered as reasonable. 
Nonetheless, we must consider the effects of picking differ- 
ent amplitudes, and we therefore consider the cases Aup - 
2.0,4.0 and 6.0 km s The motion of the piston is accompa- 
nied by luminosity variations, since the radiative flux through the 
(moving) inn er boundary is kept fix ed, as described in previous 
papers, (e.g.. iHofner & Dorfi 1997h . In order to have luminos- 
ity variations that are consistent with observations, the velocity 
amplitudes given above turn out to be a good choice. 

The inner boundary condition, i.e., the location of /?;„, was 
chosen so that all models have a similar kinetic-energy input due 
to pulsation for a given piston amplitude. This is critical in order 
to obtain a consistent pictur e of how the wind p roperties may 
vary with stellar parameters (IMattsson et al. 2008h . 



2.4. Outer boundary: Mass loss 

During the initial phase of the computation the expansion of the 
model is followed by the grid. When the outer-most point of 
the adaptive grid reaches some given distance from the inner 
boundary, the program automatically switches to a fixed outer 
boundary located at this point, allowing outflow over the outer 
boundary. There are no specific conditions posed on the location 
of the outer boundary, except that it must be at a large enough 
distance to ensure that the wind velocity at the outer boundary is 
close to the terminal velocity. At this point, the conditions 



(du\ 



0, H„ 



(4) 



are applied, where ft reflects the angular intensity distribution of 
the radiation field. 

From the models described here we obtain the density and 
the wind velocity, both as functions of radius and time. The mass 
loss rate is then a function of time given by 



M(0 = 47r/;-„,pout(f)Mout(0, 



(5) 



where pout and Mout are the density and wind velocity at Ro^t, i-e. 
the radius where the outer boundary is fixed (usually around 10'^ 
cm, or typically 20-30 stellar radii, see above). Since practically 
all momentum transfer from radiation to matter is due to the in- 
teraction with dust, we expect the wind velocity to be correlated 
with the degree of dust condensation. 
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Fig. 1. Mass-loss rates as a function of wind speeds for models 
with Aup - 4 km s ' and various values of log(C - O). 
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Fig. 2. Mass-loss rate vs. wind speed for models with diff'erent 
AMp, but the same carbon excess, log(C - O) -i- 12 = 8.80. 



2.5. Dust formation 

The models presented here include a time-dependent description 
of dust grain growth and evaporation using the moment method 
by Gail & Sedlmayr (1988) and Gauger et al. (1990). The dust 
component is described in terms of moments of the grain size 
distribution function weighted with a power j of the grain radius. 
The zeroth moment is the total number density of grains (simply 
the integral of the size distribution function over all grain sizes), 
while the third moment is proportional to the average volume 



of the grains. The moment equations are solved simultaneously 
with the RHD equations. 

In order to calculate how much momentum is transferred 
from photons to dust grains, we need to know the frequency- 
dependent opacity of these grains. This can be expressed in terms 
of the extinction efliciency, which is the ratio of the extinction 
cross section to the geometrical cross section of the grains. In the 
small-particle limit (which is used here), the dust grain opacity 
becomes a simple function of the grain radius. The dependence 
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Table 1. Definition of the grid. A denotes the grid step. 



3. Definition of the grid 
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A«p 


[Me] 




[K] 




[km s-'] 




A = 0.15 


A = 200 


A = 0.30 


A = 2.0 


0.75 


3.55 - 3.85 


2400 - 3200 


7.90-9.10 


2.0 - 6.0 


1.0 


3.70 - 4.00 


2400 - 3200 


7.90-9.10 


2.0 - 6.0 


1.5 


3.85-4.15 


2400 - 3200 


7.90-9.10 


2.0 - 6.0 


2.0 


3.85-4.15 


2400 - 3200 


7.90-9.10 


2.0 - 6.0 



of the opacity on wavelength and grain size can thus be sepa- 
rated into two independent factors, which greatly simplifies the 
calculations. However, as we shall see, this approximation may 
need to be relaxed, since the carbon grains tend to grow quite 
rapidly. The models in this paper are calculated using the refrac- 
tive index dat a of Rouleau & Martin (1991) to obtain the dust 
extinction (see I Andersen et al. 20031 : iHofner et al. 20031 for fur- 
ther discussion). The intrinsic density of the grain material used 
in the model is set to pgr = 1-85 g cm"-^, which matches the 
material in Rouleau & Martin (1991). 

We assume that dust grains can be considered to be spher- 
ical, consisting of amorphous carbo n only, s ince we are mod- 
elling carbon stars (see Andersen et al. 2003L for details about 
these assumptions). The nucleation, growth and evaporation of 
grains is assumed to proceed by reactions involving C, C2, C2H 
and C2H2. In this model of grain growth, a so called sticking co- 
efficient (sometimes referred to as the reaction efficiency factor) 
is used that enters into the net growth rate of the dust grains. 
This parameter, as, is not definitely known unless we know the 
exact sequence of chemical reactions responsible for dust forma- 
tion. However, Gail & Sedlmayr (1988) argued that the sticking 
coefficient must be of the order of unity, mainly because it is 
expected that neutral radical reactions play a major role in the 
formation of carbon grains. We set cs = 1, but note ffiat a set of 
test models with as = 0.5 give similar results. 

Dust grains in a stellar atmosphere influence both its en- 
ergy and momentum balance. For simplicity we assume com- 
plete momentum and position coupling of gas and dust, i.e. the 
momentum gained by the dust from the radiation field is di- 
rectly transferred to the gas and there is no drift between dust 
and gas. However, this strong coupling between the dust and 
the gas phase is not obvious. In a previous attempt to relax this 
phase coupling approximation, Sandin & Hofner (2003, 2004) 
found that the effects of decoupling the phases might be quite 
significant. The most striking feature is that the dust formation 
may increase significantly, but this does not necessarily increase 
ffie predicted mass loss rates for a given set of stellar param- 
eters. A coupled solution of the detailed frequency-dependent 
equations of RHD, including dust formation and drift would 
increase the computing time for each model to an extent that 
would make it very laborious to compute a grid of the kind pre- 
sented here. Thus, we chose to give frequency-dependent radia- 
tive transfer priority over relaxing the phase coupling approxi- 
mation, even if this approximation cannot hold for very low gas 
densities, i.e. at very large distances from the star or for wind 
models with very low mass loss rates. Furthermore, since trans- 
fer of internal energy between gas and dust is negligible com- 
pared to the interaction of each component with the radiative 
field (jOauger et al. 1990), we assume radiative equilibrium for 
the dust. This allows us to estimate the grain temperatures from 
the radiation temperature, as we know the dust opacities. 



Since one purpose of the grid is (among other ffiings) to provide 
information about how the mass loss rate is affected by changes 
in various stellar parameters as the star evolves on the AGB, it 
is wise to consider which values of the fundamental parame- 
ters (M, Li,, Teff) never appear in evolutionary tracks. Systematic 
studies by, e.g.. Blocker (1995) show for example that for car- 
bon stars at solar metallicity we should not expect average lu- 
minosities above 2 ■ IO^^Lq or masses above 2Mq on the AGB. 
The temperature of a carbon star lies approximately in the range 
2000 - 3500 K according to both observations (e.g. Lambert et 
al 1986, Bergeat & Chevallier 2005) and evolutionary models. 
Thus, the region in the - T^ff plane (HR-diagram) that we 
are interested in is quite small for a star of a given mass and 
initial metallicity. It is important, however, to stress that during 
its evolution an AGB star may temporarily have parameter con- 
figurations that do not 'agree' with typical observations. This is 
due to the fact that we are extremely unlikely to observe an AGB 
star as it undergoes, e.g., a thermal pulse, where the luminosity is 
very high and the temperature rather low (see Fig. 1 in Mattsson 
et al. 2007). Therefore, defining the grid according to parameter 
combinations found in evolutionary tracks will require a some- 
what larger parameter space than observations may suggest is 
necessary (see Table[T]i. 

As shown by Mattsson et al. (2008), the kinetic-energy in- 
jection by pulsation can affect the wind properties significantly. 
Placing the inner boundary at a fixed optical or geometric depth 
led to different energy injections due to different density struc- 
tures, resulting from different chemical compositions. This fact 
needs to be emphasised - theoretical work on AGB mass loss 
does not always include such "energy control" of the piston 
boundary, which may affect the results. Here we have tried 
to keep the equilibrium density/pressure at the inner boundary 
roughly the same for all models and used three different velocity 
amplitudes for the piston (Aup = 2.0, 4.0, 6.0 km s"'). 

Initially, parts of the grid were supposed to be computed us- 
ing the solar composition by Grevesse & Sauval (1998) as well 
as that of Asplund et al. (2005) , which were first thought to give 
slightly different results. However, after some scrutiny, it became 
clear that the differences were due to a poorly constrained inner 
boundary condition (see MattssoiLetal. 2008). The grid of mod- 
els presented here is computed using the Asplund et al. (2005) 
solar composition only. 



4. Results and discussion 

4.1. General trends 

Wind velocities and mass loss rates are obviously correlated 
with some stellar parameters, such as luminosity, tempera- 
ture and mass. The presence of trends with ffiese parame- 
ters can be seen in observations, but in general ffie informa- 
tion is degenerate since the different stellar parameters are 
very much coupled throughout the evolution of a star This 
is why, for e xample, one finds a simple mass-loss-period 
law (see, e.g.. Fv^ssiliadis & Wood 19931; iGroenewegen 19931 : 
iGroenewegen 19951) which of course does not contain all the in- 
formation needed for a proper treatment of mass loss in, e.g., 
stellar evolution modelling. Empirical results do not tell us di- 
rectly how the mass loss rate is affected by changes in individual 
stellar parameters. 
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Table 2. Input parameters (L*, T^ff, log(C-O), P) and the resulting avergage mass loss rate, average wind speed and the mean degree 
of dust condensation at the outer boundary for a subset of models with = IMq and Awp = 4.0km s~' . The dust-to-gas mass ratio 
Pd/Pg is calculated from /c as described in Hofner & Dorfi (1997) and the grain radius is defind as agr = ro {K\/Ko), where 
are the zeroth and first moment of the grain-size distribution, respectively. Models with log(C-0)-l-12 = 7.90 are not included, since 
none of them produced any outflow. 
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Fig. 3. Wind speed (left) and mass-loss rate (right) as functions of the acceleration parameter a for models with Aup - 4 km s~'. 
The dashed and dotted lines represent analytic models as described in Lamers & Cassinelli (1999) and Groenewegen et al. (1998), 
assuming different /M* . 



The mass-loss rate and the wind speed is, from a theoretical 
point of view, often related to the ratio of radiative to gravita- 
tional acceleration. 



r = 



AncGM^ 



(6) 



where k is the total flux-mean opacity and c is the speed of light. 
Here we define a similar quantity. 



where 

Pd _ 
Pg 



ec/c = 77 1 eo/c. 



(7) 



(8) 



in which ec - eq- sq, i.e., the carbon excess (condensible car- 
bon) and /(, is the degree of condensation, a is in principle pro- 
portional to r and can be used in its place for most purposes. The 
reason we use a instead of F is that a only contains the physi- 
cal quantities studied in this paper. In Fig. |3]we show the wind 
speed and the mass-loss rate as functions of a for all models with 
Awp = 4.0 km s ^ The wind speed is relatively well-correlated 
with a, w hich is expected from simple analytical arguments (see 
[ Lamers & Cassinel li 1999, Chap. 4 & 7). It is straight forward 
to show that, for a stationarjQ polytropic wind. 



(9) 



' The time dependent case will, at least on average, resemble the sta- 
tionary case in the sense that the wind speed is mainly due to the amount 
of dust opacity and how strong the gravitational field is. The main differ- 
ences, both in snap-shots and average structures, occur in the dust-free 
pulsating atmosphere, influencing the conditions in the dust-formation 
zone, therefore also F. 



where Wesc is the escape speed, Cs is the sound speed, is the 
condensation radius and y is the polytropic index in a polytropic 
equation of state. For the mass-loss rate, the correlation with F 
(or a) is weaker, which is due to variations of the gas-density 
profiles of the wind regions between different models. However, 
the mass-loss rate can (in the stationary case) be expressed an- 
alytically as (see Chap. 7 in Lamers & Cassinelh 1999, and the 
Appendix in Groenewegen et al. 1998) 



M OCT^ 1 1 - - 

Mout - u{Rc) \ F 



(10) 



where t„ is the flux-mean optical depth far out in the wind, and 
all other symbols are as previously defined. By analysing the 
time development of several wind models, it becomes obvious 
that the density profiles are not directly a consequence of the 
momentum transfer, as the wind speed profiles are, but of a com- 
bination of momentum transfer and pulsation dynamics. Thus, 
one should not expect a strong coiTelation between mass-loss 
rate and a for the non-stationary case. 

Previous efforts to model the winds of carbon stars are 
in many cases not directly comparable to the present results. 
For instance, stationary wind models with detailed dust for- 
mation should not be compared with DMAs, nor should mod- 
els including drift but having grey radiative transfer be com- 
pared with models not including drift but featuring frequency- 
dependent radiative transfer instead (as in the present work). In 
comparison with existing theoret ical work similar to the present 
(Winters et al. 2000; Wachter et al. 20021) . our mass-loss rates 
show over-all trends that are comparable, although individual 
combinations of stellar parameters in some cases lead to quite 
different results c ompared to their mo dels. More recent work by 
the Berlin group ( Wachter et al. 20081) has been focusing on sub- 
solar metallicities and is not adequate for comparison. The wind 
speeds that we obtain for the most carbon rich m odels are higher 
than the highest observed (cf. lSchoier et al. 200 ID . but consistent 
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Fig. 4. Wind speed (left) and mass-loss rate (right) as functions of the dust-to-gas ratio for models with Awp = 4 km s ' . Grey dots 
shows modelling results by Winters et al. (2000) and the dashed line marks the division between A- and B-type models according 
to Winters et al. (2000). 



with the results by Winters (2000) and Wachter et al. (2002) for 
similar dust-to-gas ratios (see Fig. |4]i. 

4.2. Wind velocity 

The wind velocity depends strongly on the carbon excess sc, 
which appears as a steep slope in Fig.|6l This is expected, since 
the acceleration of the gas depends on the dust opacity, which 
in turn is proportional to the dust-to-gas ratio. The strong cor- 
relation between dust abundance and wind velocity (see Fig.|4]i 
demonstrates how the wind speed is a direct consequence of mo- 
mentum transfer from radiation to gas via dust grains. In addi- 
tion, that the abrupt decline in the wind speeds at high tempera- 
tures (see Fig. |6]) is mostly due to inefficient dust formation. At 
high Teff the condensation radius is moved out to a distance 
too far away from the star for any significant amount of material 
to reach this radius. Consequently, almost no dust is formed and 
therefore no wind either 

4.3. Mass-loss rate 

The classical Reimers (1975) formula for the mass loss rate of 
red giants, 

^3/2 

Mr ^ 1.3410-^7; * , , (11) 

where we have used - Ano-R^T^^ to eliminate can 
be obtained by relating the kinetic energy of the outflow to 
the radiative energy flux and the gravitational potential en- 
ergy. The Reimers law was derived strictly for red giants and 
cannot be applied to AGB stars without significant modifica- 
tions (^udritzki & Reimers 1978). Fitting a simple function to 
a set of wind models reveals a stronger dependence on lu- 
minosity and effective temperature for carbon stars (see, e.g.. 



lArndt et al. 19971 IWachter et al. 20021) . In Fig. |9] one may note 
a weak correlation between the mass-loss rates obtained from 
our detailed modelling and the Reimers formula, but since the 
Reimers formula (and commonly used modified versions of it, 
e.g., the formula suggested by Blocker 1995) does not include 
any dependence on the free carbon abundance, it cannot easily 
be modified so that it would reproduce our results. Using a mod- 
ified Reimers law to account for the mass-loss during the carbon 
star phase in stellar evolution models is therefore not to be rec- 
ommended. 

In general, the mass loss rate follows a slightly different pat- 
tern than the wind velocity. First of all there is an obvious anti- 
correlation between effective temperature Teff and the mass loss 
rate (see Fig.|6]l, which partly can be interpreted as an effect of 
the Teff - \og{g) correlation. But this anti-correlation is not di- 
rectly due to the effect of surface gravity, log(g), on the wind dy- 
namics. With increasing log(g), the star becomes more compact 
and the outer atmospheric layers become thinner, which means 
that there will be less material available in the dust-formation 
zone around the condensation radius. With lower density where 
the dust forms, the mass loss rate may drop, since there will be 
less gas available in the wind forming region. 

The parameter configurations we have used here give mass- 
loss rates mostly within the interval 10"^ - lO^^M© yr"', which 
is in agree ment with mass los s rates derived from observations 
(see, e.g., ISchoier et al. 200 iL and Fig. [T] in this paper), apart 
from the lack of models with mass loss rates below lO^^Moyr"' 
and the correlation between mass-loss rate and wind speed ob- 
tained from observations. This correlation could be an effect of 
the coupling between stellar parameters as the stars evolve on 
the AGB. But a detailed discussion about how to compare the 
present model grid with observations so that biases from sparse 
grid sampling and observational selection effects are reduced to 
a minimum, is beyond the scope of this paper However, we may 
note that stars will stay much longer in evolutionary states that 
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correspond to low mass-loss rates than those that correspond the 
superwind phase (typically the tip of the AGB) during which the 
mass-loss rate is very high and the evolution is rapid. Hence, one 
is more likely to find a star with a low mass loss rate. The low 
mass-loss rates may correspond to the transition regions where 
the mass-loss thresholds appear, and these transitions from no 
(or very inefficient) mass loss to strong, dust driven winds are 
perhaps not adequately resolved, the way the grid is defined. We 
have explored some of these transition regions in more detail 
by taking smaller steps in effective temperature and luminosity, 
which showed that the wind appears to be "switched on" at more 
or less a specific parameter value (see Fig. [8]). Whether this steep 
threshold is a real property of dust driven winds or a model arte- 
fact is not clear at present, although mass-loss thresholds should 
exist for physical reasons (see Sect. 14.41 ). It is possible that a 
smooth transition is more realistic and that the observed stars 
with low mass-loss rates are sitting in these transition regions, 
which may correspond to very restricted parameter intervals. We 
will return to this important issue in a forthcoming paper. 

The energy injection by pulsations, or more precisely, the 
kinetic energy density at the inner boundary of the model atmo- 
sphere, also plays an important role in the wind-formation pro- 
cess. A correlation between this quantity and the mass loss rate 
obtained for specific sets of stellar parameters was pointed out by 
Mattsson et al. (2008) and stems from the fact that the pulsations 
are needed to levitate the atmosphere such that dust formation 
can take place. The efficiency of dust formation is therefore more 
or less affected by the pulsation energy. Hence, for AGB stars, it 
is necessary to consider the kinetic energy from the pulsations as 
a part of the wind-driving mechanism, although it is the momen- 
tum transfer from radiation to dust that maintains the outflow. 
We have tried different piston-velocity amplitudes and, as ex- 
pected, stronger puls ations seem to favour wind formation (see 
iMattsson et al. 20081 for further details), although such effects 
of AMp are only important in the critical wind regime (near the 
mass-loss threshold). The most obvious effect of changing AMp is 
that reducing the strength of the pulsations means that a smaller 
region in stellar parameter space will correspond to stable out- 
flows and vice versa, i.e., changing Awp can shift the loca tions 
of the mass-loss thresholds (see Fig. 2 in lMattsson et al. 2007bl) . 
but the typical mass-loss rate as well as the range of mass-loss 
rates in the wind-forming cases are not affected significantly (see 
Fig. 13. 

Comparing the mass-loss rates we obtain with the "single- 
scattering limit" mass loss, i.e.. 



Mo = 



1 



(12) 



where c is the speed of light, we find that a fair number of models 
where the mass-loss rate is higher than Mq, which can be seen 
in Fig.|5] The quantity Mq is, however, not a proper upper limit 
for "single scattering" due to the assumptions made in deriving 
it (see, e.g., Lamers & Cassinelli 1999). The optical depth of the 
wind may very well be larger than unity around the flux peak in 
some cases, even if the total optical depth is not. Furthermore, 
it has been shown by Gail & Sedlmayr (1986) that heavily ob- 
scured objects must have mean optical depths that exceed unity, 
so the result in Fig.|5]is hardly surprising. Furthermore, the time- 
dependent luminosity L may not be represented properly by the 
luminosity L* of the static start model used in Eq. ( fT2b and in 
Fig El since L is phase dependent in the dynamical models and 
wind acceleration may occur within a limited range of phases in 
a given model, thus changing the "effective luminosity. 
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Fig. 5. Modelled mass-loss rates vs. the "single-scattering limit" 
(see Eq. [T2l l. The dashed line marks the case of a one-to-one 
correspondance. 

4.3.1. Dust-to-gas ratios 

The dust opacity scales with the relative dust abundance Pd/Pg- 
Thus, we should expect the wind velocity to behave in a way 
similar to Pd/Pg, which is actually the case if one compares Fig. [6] 
with Fig.[7l According to ob servations (e.g.. lSchoier et al. 200 iT 
IBergeat & Chevalher 2005b . the pd/pg-ratio is nearly indepen- 
dent of Teff. This fairly consistent with the weak Teff-dependence 
that we obtain. However, the clearly increasing trend with the 
abundanc e of condensible carbon eq is not seen in observations 
(see, e.g.. lBergeat & Chevallier 2005h although the existence of 
such a trend should be rather obvious since pd o^eq. 

Because the stellar luminosity and mass is kept constant in 
Fig. |2] the plotted dust-to-gas ratio will also tell us how the ac- 
celeration ratio a varies with Teff and the abundance of free car- 
bon, ec- Below a certain T^s, the dust-to-gas ratio pd/Pg depends 
mainly on ec, and we may thus conclude that a is to a large ex- 
tent set by Ec for given a mass and luminosity. 

4.3.2. Average grain size 

We return now to the grain-size issue and whether the small- 
particle limit actually applies to the dusty atmospheres of carbon 
stars. The characteristic grain radius can be computed as 



^^1 



(13) 



where ro is the monomer radius and Ao, K\ are the zeroth and 
first moment of the grain size distribution function, respectively. 
Using this definition of grain radius, we obtain a typical grain 
radius of Ogi ~ lO"*" cm in the set of models with M* = IMq, 
AMp - 4.0 km s that produce winds. The small particle ap- 
proximation is assumed to be applicable as long as the true radi- 
ation pressure efficiency factor 2ip does not deviate from that of 
the approximation by more than 10% at 1/zm. Thus, the small- 
particle limit requires that < 10^^ '^ cm, i.e, not larger than a 
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Fig. 6. Mean mass-loss rate (upper panels), mean dust-loss rate (middle panels) and wind speed (lower panels) as functions of 
effective temperature and the abundance of condensible carbon for models with - I.OMq and Aup - 4 km s 



few times 10"*' cm (see Fig.fTTTi. which indicates that the grains 
forming in our models are on average a little too large for 2rp 
to be well-approximated by the small-particle limit. The wind 
speed is anti-correlated with dust-grain size (see Fig.fTOli. which 
we interpret as a simple consequence of the fact that for a slow 
wind the grains stay longer in the dust-formation zone and thus 
they have time to grow bigger. 

Since 2,p (which describes the net effect of momentum 
transfer by absorption and scattering for dust) will level-out to 
a constant above a certain grain radius, i.e., 2rp is only propor- 
tional to the grain size for small grains, our model (which is 
based on the small particle approximation of Q^p, irrespective 



of grain size) will tend overestimate the dust opacity for really 
large particles. However, in the transition region, where Q^p goes 
from being proportional to the grain size to being almost inde- 
pendent of the grain size, the small-particle approximation will 
instead underestimate the extinction by as much as a factor of 
five (again, see Fig. [TTT l. Since there is a competition between 
nucleation and grain growth, we expect this underestimation to 
be important only in critical cases with slow winds. The majority 
of the grid models, however, show grain sizes that correspond to 
the transition region where is likely underestimated. Exactly 
how this may affect the mass-loss rates and other results in this 
paper is difficult to estimate and needs further investigation. 
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Fig. 7. Mean dust-to-gas ratio (upper panels) and degree of dust condensation (lower panels) as functions of effective temperature 
and the abundance of condensable carbon for models with M^, - 1 .OMq and Aur, - 4 km s ' . 
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4.4. Mass-loss thresholds 

Mattsson et al. (2007b) demonstrated the necessity of consider- 
ing thresholds in mass-loss prescriptions for dust-driven winds. 
As we mentioned above, a mass-loss threshold appears, as one 
would expect, and we find that for ec < 8.20 and/or T^ff > 3200 
K a dust driven does not seem to form regardless of the com- 
bination of other stellar parameters (see Fig. |6] and Fig. |7] as 
well as the additional online material). Since the surface gravity 
scales with T^^fi and dust formation is a process highly sensitive 
to temperature, it is expected that the mass-loss rate is highly 
dependent on Teff too. But, as shown above, there is a rather 
strong dependence on condensible carbon for both the wind ve- 
locity and the mass-loss rate, which is quite interesting in com- 
parison with previous studies of this kind. Arndt et al. (1997) as 
well as Wachter et al. (2002) find a rather weak dependence on 
C/O (which is their choice of parameter). Their results stand in 
quite sharp contrast to the results presented here. However, our 
findings here are, qualitatively speaking, hardly a new discov- 
ery. Hofner and Dorfi (1997) and Winters (2000) have already 
pointed out the C/O-dependence, especially in the critical wind 
regime, although this has not been widely recognised. 

The mass-loss thresholds originate from the simple fact that 
there exists a critical acceleration ratio aa-it and that there con- 
sequently also exist critical values for sq, Teff and L* for which 
a = aciit- The principles behind this may easily be understood 
from basic physics, but none of the existing parametric formu- 
lae for dust-driven mass loss (e.g. Arndt et al. 1997, Wachter 
et al. 2002) contain any thresholds. These transition regions are 
perhaps not adequately covered due to the grid spacing we have 
chosen and the sudden onset of an outflow discussed earlier may 
partly be due to model assumptions. But, the existence of mass- 
loss thresholds as such is a fact, and should definitely be included 
in a mass-loss prescription for stellar evolution modelling, even 
if it is unlikely that one can find a sufficiently simple functional 



form that will cover all aspects of how mass loss depends on 
stellar parameters. 

The threshold due to the amount of carbon excess implies 
that Galactic carbon stars with low observed C/O ratios should 
typically show no or very little mass loss. Some support for this 
idea comes from the work by Lambert et al. (1986), where about 
50% of their sample of carbon stars showed C/O < 1 . 1 (all stars 
in their sample have roughly solar oxygen abundances), and a 
fair number of those stars do not have any detected outflows. A 
puzzeling fact, howe ver, is that many s tars with C/O ~ 1 do have 
rather strong winds dRamstedt et al. 20 06), observations which 
cannot be properly explained at present. 

4.5. No simple formula! 

Consider two stars with the same mass, luminosity and effective 
temperature. Increasing the amount of condensible carbon ec by 
a factor of 5 may increase the mass loss rate by almost an order 
of magnitude (see Fig.|6]l. This may not appear as a huge prob- 
lem in a log-log-plot, as in Fig.|9]where the grid models are com- 
pared to a scaled Reimers (1975) law, even if ec is not a param- 
eter in Reimers' formula. But note that there is a 0.5 dex scatter 
in Fig. |9] (dotted lines) which actually makes a considerable dif- 
ference for, e.g., stellar evolution modelling, where ec changes 
during the carbon-star phase, causing a different mass-loss evo- 
lution and the occurrence of a superwind. A direct fit to the 
whole set of data, including both ec and Awp as fitting parameters 
does not really fix the problem, since statistically preferred cor- 
relations (formulae obtained through, e.g., least-square fitting or 
;if^-minimisation) on ly show a rather weak dependence on t hese 
two parameters (see lArndt et al. 19971: [Wachter et al. 2002 h. On 
the other hand, if one keeps all other parameters fixed and then 
varies ec or Awp, the mass-loss rate is definitely strongly affected 
and it therefore seems dangerous to conclude that these two pa- 
rameters can be omitted in the prescription as Arndt et al. (1997) 
and Wachter et al. (2002) suggested. 

For the reasons given above, we do not recommend that for- 
mulae derived from the data given in this paper are used to re- 
place the actual grid. Instead, using the tables (on-line material) 
to create a look-up matrix for the mass losfl (combined with 
multi-dimensional interpolation) in stellar evolution modelling 
will guarantee that detailed features (such as the increase of the 
mass loss rate as the amount of condensible carbon increases) 
will be properly included, while simplified formulae will not. As 
an example, Mattsson et al. (2007) showed that, as the amount of 
condensible carbon changes rapidly (relatively speaking, often 
more rapidly than luminosity, effective temperature etc.) during 
the TP-AGB phase, simple prescriptions of mass loss may there- 
fore not be correct, since that parameter (condensible carbon) is 
usually not included. Note, however, that the mass-loss prescrip- 
tion we provide here best describes the superwind of carbon stars 
(at the tip of the AGB) and winds corresponding to the early/mid 
AGB (where mass-loss rates are much lower) may not always be 
accurately covered (see Sect. 14.41 and below). 

4.6. Uncertainties 

Generally, numerical models are subject to various uncertain- 
ties due to simplifying assumptions, numerical limitations and 
the fact that some relevant aspects of the physical reality may 
not be included in the model. Several of the models presented 



- A FORTRAN routine for this purpose is available on-line at 
http://coolstars.astro.uu.se| 
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Fig. 10. Wind speed (left) and mass-loss rate (right) versus the characteristic grain radii for models with Auj, = 4 km s ' and 
Mi, - IMq. Note the absence of any correlatiopn with mass-loss rate. 



here show intermittent, chaotic behaviour that may, or may not, 
reflect non-linearities in the physics of real dynamic stellar at- 
mospheres. The existence of such non-linear features can affect 
the time averages of mass-loss rates, wind speeds and other wind 
properties that we have derived here. By necessity, the time se- 
ries from which these averages are computed must be limited 
and the time series that we obtain for a specific set of stellar pa- 
rameters over the 100-400 pulsation periods the models cover 
may not always be representative. For example, the "deviating" 
mass-loss rate seen in Fig.|6]for log{Li, / Lq) - 4.00, T = 3000 
K and log(C-O) + 12 - 9.10 is probably an effect of dynamical 
non-linearities, which turn out to be favourable for dust forma- 
tion. But it is hard to know for certain without a very long time 
series, and the models cannot be evolved for more than a few 
hundred periods, since it would lead to significant mass deple- 
tion in the modelled region. For the more "well-behaved" cases, 
the relative error in the time average due to time variations of the 
mass-loss rate, wind speed etc., is typically no bigger than 10% 
(usually much less). 

Another source of uncertainty (or bias) is our choices of stel- 
lar parameters, i.e., the definition of the grid. As previously dis- 
cussed, the threshold regions may not be adequately covered, 
which can affect the interpolations in the mass-loss prescrip- 
tion described above. As shown in Fig. [8] linear interpolation is 
probably sufficiently accurate in most cases. But the very steep 
thresholds we find when the grid spacing is decreased may not 
be reflecting the true nature of these critical regions. It cannot be 
excluded that one possible explanation for the steepness might 
be the resolution of the radiative transfer, which is relatively low 
(64 frequency points) in the present study. In critical cases, a 
higher resolution may slightly affect the momentum transfer effi- 
ciency (from radiation to dust). We have recently begun to study 
the effects of frequency resolution (i.e., computing models with 
high-resolution radiative transfer), but it is yet too early to draw 
general conclusions. It is also possible that the effects of dust 
grain sizes can affect the radiative transfer significantly in mod- 



els near the thresholds. These effects are currently being investi- 
gated as well, and will be discussed in a forthcoming paper. 

Finally, there are two, possibly important, ingredients that 
are not included the present DMAs: dynamical gas/dust decou- 
pling (i.e., drift) and the effects of the size distribution of dust 
grains on the momentum transfer effic iency. Grain size appears 
to be a critical issue for M-type stars (Hofner 2008) and it may 
be of some importance also for the winds of carbon stars. It is 
kno\ yn that drift does have a r i effect on the wind prop erties (see, 
e.g.. ISandin & Hofner 2001: ISandin & Hofner 2004 . but it has 
not been studied for the frequency-dependent case together with 
a detailed description of dust condensation in time-dependent 
wind models. 



5. Summary and conclusions 

The present grid of wind models has shown that mass-loss rates 
and wind speeds are clearly affected by the choice of stellar tem- 
perature, mass, luminosity and the abundance of available car- 
bon. Furthermore, there are also inevitable mass-loss thresholds 
below which a dust-driven wind cannot form, in certain parts 
of the parameter space. Contrary to some previous studies, we 
find a strong dependence on the abundance of free carbon, which 
turns out to be a critical parameter. In a stellar evolution context, 
we expect that this dependence on the carbon excess will natu- 
rally lead to the development of a superwind after a few dredge- 
up events. Hence, the AGB would be terminated as soon as the 
atmosphere has become sufficiently carbon rich, and the number 
of thermal pulses will be limited by a self -regulating mechanism: 
each thermal pulse increases the carbon excess, which in turn in- 
creases the mass loss. 

The set of models with Aup = 4.0 km s"' show mass-loss 
rates that are typically a few times 10^^ Mq yr"' and only one 
model has a really low mass-loss rate M < 10"^. With Aup - 2.0 
km s"' and Auj, = 6.0 km s"' the region in the stellar parameter 
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amorphous carbon dust evalu- 
ated at A = lyum. The vertical, 
thin, dashed line shows where 
the small particle approxima- 
tion breaks down. 



space where winds occur is somewhat smaller or larger, respec- 
tively. The typical mass-loss rates are otherwise similar to that 
of the Awp = 4.0 km s"' subset. Models with a steady mass 
loss rate M < 10"^ appear to be difficult to produce regardless 
of the combination of stellar parameters and piston amplitude. 
The difficulty of producing really low mass-loss rates may rep- 
resent a shortcoming of the present wind model. It is yet unclear 
whether this is due to, e.g., the low-resolution radiative transfer 
of the models, or the phase coupling between gas and dust (i.e., 
absence of drift). A smaller number of test models with high 
resolution (440 frequency points instead of 64) has been com- 
puted, but no dramatic differences were found. Further study of 
this matter is needed, however. 

Winds of carbon stars are known to be relatively fast, com- 
pared to other long-period variables. The wind speeds that we 
obtain lie in the range 1 - 60km s"', where the fastest winds are 
found in the most carbon rich cases, i.e., carbon stars on the tip 
of the AGB. The lack of observed carbon stars with very high 
wind speeds (mqui > 40km s"') may be a selection effect, since 
the chance of finding a carbon star at the very final stage of its 
evolution must be quite small. The very low wind speeds seen in 
a few models (uout < 5km s"'), but not in observed carbon stars, 
is an interesting prediction of the models. Such slow outflows 
may exist, since a wind can, in theory, be maintained as long as 
the flow surpasses the escape speed at some point in the velocity 
profile. 

The mean degree of dust condensation in the wind-producing 
models varies between ( fc) ^ 0.1 and ( fc) ~ 0.7, while the typ- 
ical value is around 0.3. In the light of the recent models of M- 
type stars by Hofner (2008) , we derived the characteristic grain 
sizes for a subset of models. The characteristic grain radius in 
the the models that produce winds is Ogr ~ 10"^' cm, which is a 
somewhat too large for the small-particle limit to be a good ap- 
proximation. The dust-grain size shows an anti-correlation with 



the wind speed, while the mass-loss rate seems totally uncorre- 
cted. Large grains tend to form in slow winds, since they stay 
longer in the dust-formation zone and therefore have extra time 
to grow. However, the anti-correlation seen in Fig.[TO]may not be 
entirely correct due to the fact that we compute the dust extinc- 
tion in the small-particle limit, which is not fully consistent with 
the grain sizes we obtain from the model of dust formation. The 
grain size is also anti-correlated with the abundance of condensi- 
ble carbon, which is likely an effect of the competition between 
nucleation and grain growth. 

Despite some shortcomings in reproducing low mass-loss 
rates, the present grid is still a significant step forward. This 
is the first systematic study of mass loss and other wind prop- 
erties as functions of stellar parameters that is based on a dy- 
namic atmosphere model including both frequency-dependent 
radiative transfer and detailed, time-dependent dust formation, 
and is probably the most advanced theoretical constraint on the 
super-wind of carbon stars at present. Our results have shown 
that strong, dust driven winds cannot develop under arbitrary 
physical circumstances, i.e., in certain regions of (stellar) param- 
eter space, carbon stars do not experience any significant dust- 
driven mass loss. 

The thresholds that appear when, e.g., the effective temper- 
ature becomes too high or the abundance of condensible carbon 
too low, must be taken into account when constructing a mass- 
loss prescription to be used with other types of stellar models, 
such as stellar evolution and nucleosynthesis. Simple mass-loss 
formulae cannot reproduce all aspects of these thresholds and 
we have therefore not derived any such formula. Instead, we pro- 
vide an easy-to-use FORTRAN-code that finds the correspond- 
ing mass-loss rate from interpolated data cubes. The code will 
soon be made available on-line. 
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Table 3. Input parameters (L*, T^g, C/O, P) and the resulting mean mass loss rate, mean velocity at the outer boundary and mean 
degree of dust condensation at the outer boundary for a subset of models with = O.TSMg and Aup = 2.0km s~' . The dust-to-gas 
mass ratio Pdust/Pgas is calculated from /c as described in Hofner & Dorfi (1997). Where "***" is given instead of a number, no 
meaningful wind properties could be derived due to numerical problems. 



log(L*) 


^eff 


log(C-0)+12 


P 




("out) 


(fc> 


(Pd/Pg> 




[K] 




[days] 


[Moyr-i] 


[kms-'] 






3.55 


2400 


8.20 


221 










3.55 


2400 


8.50 


221 


3.95E-07 


4.46E+00 


2.38E-01 


6.45E-04 


3.55 


2400 


8.80 


221 


1.62E-06 


1.86E+01 


2.89E-01 


1.56E-03 


3.55 


2400 


9.10 


221 


2.53E-06 


3.70E+01 


4.51E-01 


4.87E-03 


3.70 


2400 


8.20 


295 










3.70 


2400 


8.50 


295 


1.46E-06 


9.01E+00 


2.21E-01 


5.99E-04 


3.70 


2400 


8.80 


295 


3.62E-06 


2.13E+01 


3.11E-01 


1.68E-03 


3.70 


2400 


9.10 


295 


4.62E-06 


3.93E+01 


4.80E-01 


5.18E-03 


3.85 


2400 


8.20 


393 


1.06E-06 


1.68E+00 


2.28E-01 


3.10E-04 


3.85 


2400 


8.50 


393 


4.59E-06 


1.29E+01 


2.44E-01 


6.61E-04 


3.85 


2400 


8.80 


393 


6.81E-06 


2.46E+01 


3.37E-01 


1.82E-03 


3.85 


2400 


9.10 


393 


8.09E-06 


4.31E+01 


5.34E-01 


5.76E-03 


3.55 


2600 


8.20 


221 










3.55 


2600 


8.50 


221 


3.13E-07 


6.79E+00 


2.18E-01 


5.91E-04 


3.55 


2600 


8.80 


221 


9.67E-07 


2.28E+01 


2.96E-01 


1.60E-03 


3.55 


2600 


9.10 


221 


1.46E-06 


4.10E+01 


3.88E-01 


4.19E-03 


3.70 


2600 


8.20 


295 










3.70 


2600 


8.50 


295 


7.44E-07 


6.86E+00 


1.69E-01 


4.58E-04 


3.70 


2600 


8.80 


295 


2.34E-06 


2.35E+01 


2.70E-01 


1.46E-03 


3.70 


2600 


9.10 


295 


2.83E-06 


4.24E+01 


3.95E-01 


4.26E-03 


3.85 


2600 


8.20 


393 


4.18E-07 


2.07E+00 


2.01E-01 


2.73E-04 


3.85 


2600 


8.50 


393 


2.03E-06 


l.OlE+01 


1.56E-01 


4.23E-04 


3.85 


2600 


8.80 


393 


*** 


*** 


*** 


*** 


3.85 


2600 


9.10 


393 


5.90E-06 


4.31E+01 


4.26E-01 


4.60E-03 


3.55 


2800 


8.20 


221 










3.55 


2800 


8.50 


221 










3.55 


2800 


8.80 


221 










3.55 


2800 


9.10 


221 


1.33E-06 


3.80E+01 


3.82E-01 


4.12E-03 


3.70 


2800 


8.20 


295 










3.70 


2800 


8.50 


295 










3.70 


2800 


8.80 


295 


8.62E-07 


2.57E+01 


2.24E-01 


1.21E-03 


3.70 


2800 


9.10 


295 


2.72E-06 


4.33E+01 


4.39E-01 


4.74E-03 


3.85 


2800 


8.20 


393 










3.85 


2800 


8.50 


393 


1.36E-06 


1.13E+01 


1.38E-01 


3.74E-04 


3.85 


2800 


8.80 


393 


2.59E-06 


2.69E+01 


2.48E-01 


1.34E-03 


3.85 


2800 


9.10 


393 


3.73E-06 


4.70E+01 


3.84E-01 


4.14E-03 


3.55 


3000 


8.20 


221 










3.55 


3000 


8.50 


221 










3.55 


3000 


8.80 


221 










3.55 


3000 


9.10 


221 










3.70 


3000 


8.20 


295 










3.70 


3000 


8.50 


295 










3.70 


3000 


8.80 


295 










3.70 


3000 


9.10 


295 


1.23E-06 


3.80E+01 


2.87E-01 


3.10E-03 


3.85 


3000 


8.20 


393 










3.85 


3000 


8.50 


393 










3.85 


3000 


8.80 


393 


1.14E-06 


2.71E+01 


2.21E-01 


1.20E-03 


3.85 


3000 


9.10 


393 


3.44E-06 


4.59E+01 


3.92E-01 


4.23E-03 


3.55 


3200 


8.20 


221 










3.55 


3200 


8.50 


221 










3.55 


3200 


8.80 


221 










3.55 


3200 


9.10 


221 










3.70 


3200 


8.20 


295 










3.70 


3200 


8.50 


295 










3.70 


3200 


8.80 


295 










3.70 


3200 


9.10 


295 


2.47E-08 


1.26E+01 


3.64E-01 


3.93E-03 


3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 










3.85 


3200 


9.10 


393 
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Table 4. Same as[3] but for Awp = 4.0km s"' . 



log(L*) 




log(C-0)+12 


P 




<"oiit> 


<fc> 


iPd/Pg) 


[Lo] 


[K] 




[days] 


[Moyr-'] 


[km s-'] 






3.55 


2400 


8.20 


221 


- 


- 


- 


- 


3.55 


2400 


8.50 


221 


1.55E-06 


5.80E+000 


2.89E-01 


7.83E-04 


3.55 


2400 


8.80 


221 


3.05E-06 


1.98E+001 


3.80E-01 


2.06E-03 


3.55 


2400 


9.10 


221 


3.39E-06 


3.66E+001 


5.10E-01 


5.50E-03 


3.70 


2400 


8.20 


295 


- 


- 


- 


- 


3.70 


2400 


8.50 


295 


2.77E-06 


8.40E+000 


2.32E-01 


6.29E-04 


3.70 


2400 


8.80 


295 


4.83E-06 


2.11E+001 


3.32E-01 


1.80E-03 


3.70 


2400 


9.10 


295 


5.84E-06 


3.84E+001 


5.11E-01 


5.51E-03 


3.85 


2400 


8.20 


393 


2.23E-06 


3.32E-01 


3.07E-01 


4.17E-04 


3.85 


2400 


8.50 


393 


7.72E-06 


1.18E+001 


2.46E-01 


6.67E-04 


3.85 


2400 


8.80 


393 


9.07E-06 


2.36E+001 


3.49E-01 


1.89E-03 


3.85 


2400 


9.10 


393 


1.05E-05 


4.16E+001 


5.69E-01 


6.14E-03 


3.55 


2600 


8.20 


221 


- 


- 


- 


- 


3.55 


2600 


8.50 


221 


4.35E-07 


5.03E+000 


2.19E-01 


5.94E-04 


3.55 


2600 


8.80 


221 


1.26E-06 


2.06E+001 


2.74E-01 


1.48E-03 


3.55 


2600 


9.10 


221 


1.97E-06 


3.89E+001 


4.13E-01 


4.46E-03 


3.70 


2600 


8.20 


295 


- 


- 


- 


- 


3.70 


2600 


8.50 


295 


1.33E-06 


7.44E+000 


1.86E-01 


5.04E-04 


3.70 


2600 


8.80 


295 


2.84E-06 


2.21E+001 


2.65E-01 


1.43E-03 


3.70 


2600 


9.10 


295 


4.13E-06 


3.96E+001 


4.36E-01 


4.70E-03 


3.85 


2600 


8.20 


393 


l.OlE-06 


6.21E-01 


2.19E-01 


2.98E-04 


3.85 


2600 


8.50 


393 


4.11E-06 


1.13E+001 


1.96E-01 


5.31E-04 


3.85 


2600 


8.80 


393 


3.99E-06 


2.82E+001 


3.43E-01 


1.86E-03 


3.85 


2600 


9.10 


393 


5.66E-06 


4.65E+001 


5.25E-01 


5.67E-03 


3.55 


2800 


8.20 


221 


- 


- 


- 


- 


3.55 


2800 


8.50 


221 


- 


- 


- 


- 


3.55 


2800 


8.80 


221 


7.18E-07 


2.56E+001 


3.19E-01 


1.73E-03 


3.55 


2800 


9.10 


221 


1.87E-06 


4.12E+001 


4.82E-01 


5.20E-03 


3.70 


2800 


8.20 


295 


- 


- 


- 


- 


3.70 


2800 


8.50 


295 


8.69E-07 


8.85E+000 


1.69E-01 


4.58E-04 


3.70 


2800 


8.80 


295 


1.77E-06 


2.58E+001 


2.83E-01 


1.53E-03 


3.70 


2800 


9.10 


295 


2.84E-06 


4.45E+001 


4.56E-01 


4.92E-03 


3.85 


2800 


8.20 


393 


- 


- 


- 


- 


3.85 


2800 


8.50 


393 


2.47E-06 


1.03E+001 


1.44E-01 


3.90E-04 


3.85 


2800 


8.80 


393 


3.57E-06 


2.69E+001 


2.25E-01 


1.22E-03 


3.85 


2800 


9.10 


393 


4.77E-06 


4.68E+001 


5.14E-01 


5.55E-03 


3.55 


3000 


8.20 


221 


- 


- 


- 


- 


3.55 


3000 


8.50 


221 


- 


- 


- 


- 


3.55 


3000 


8.80 


221 


- 


- 


- 


- 


3.55 


3000 


9.10 


221 


1.06E-06 


3.71E+001 


3.81E-01 


4.11E-03 


3.70 


3000 


8.20 


295 


- 


- 


- 


- 


3.70 


3000 


8.50 


295 


- 


- 


- 


- 


3.70 


3000 


8.80 


295 


9.35E-07 


2.61E+001 


2.81E-01 


1.52E-03 


3.70 


3000 


9.10 


295 


2.67E-06 


4.31E+001 


4.65E-01 


5.02E-03 


3.85 


3000 


8.20 


393 


- 


- 


- 


- 


3.85 


3000 


8.50 


393 


1.07E-06 


1.21E+001 


1.35E-01 


3.66E-04 


3.85 


3000 


8.80 


393 


2.34E-06 


2.72E+001 


2.14E-01 


1.16E-03 


3.85 


3000 


9.10 


393 


3.85E-06 


5.03E+001 


4.78E-01 


5.16E-03 


3.55 


3200 


8.20 


221 


- 


- 


- 


- 


3.55 


3200 


8.50 


221 


- 


- 


- 


- 


3.55 


3200 


8.80 


221 


- 


- 


- 


- 


3.55 


3200 


9.10 


221 


- 


- 


- 


- 


3.70 


3200 


8.20 


295 










3.70 


3200 


8.50 


295 










3.70 


3200 


8.80 


295 


5.32E-07 


2.13E+001 


1.27E-01 


6.87E-04 


3.70 


3200 


9.10 


295 


1.56E-06 


4.64E+001 


3.11E-01 


3.36E-03 


3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 


1.05E-06 


2.44E+001 


2.10E-01 


1.14E-03 


3.85 


3200 


9.10 


393 


3.92E-06 


4.70E+001 


5.00E-01 


5.40E-03 
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Table 5. Same as[3] but for Awp = 6.0km s"' . 



log(L*) 




log(C-0)+12 


P 




<"oiit> 


<fc> 


iPd/Pg) 


[Lo] 


[K] 






[days] 


[Moyr-'] 


[km s-'] 






3.55 


2400 








221 


- 


- 


- 


- 


3.55 


2400 








221 


1.35E-06 


3.08E+000 


2.58E-01 


6.99E-04 


3.55 


2400 


Q 

o 


. ou 


221 


2.86E-06 


1.71E+001 


3.07E-01 


1.66E-03 


3.55 


2400 


q 

y 


1 fl 


221 


4.11E-06 


3.53E+001 


5.27E-01 


5.69E-03 


3.70 


2400 


Q 
O 


9n 


295 


2.57E-06 


9.87E-01 


4.77E-01 


6.48E-04 


3.70 


2400 


Q 
O 


Sf) 

. JU 


295 


3.99E-06 


6.13E+000 


2.22E-01 


6.02E-04 


3.70 


2400 


a 
o 


. ou 


295 


6.54E-06 


2.06E+001 


3.50E-01 


1.89E-03 


3.70 


2400 


q 

y 


1 fl 


295 


6.99E-06 


3.79E+001 


5.50E-01 


5.93E-03 


3.85 


2400 


a 
o 




393 


- 


- 


- 


- 


3.85 


2400 


a 
o 


SO 


393 


1.05E-05 


9.70E+000 


2.11E-01 


5.72E-04 


3.85 


2400 


a 
o 


80 


393 


1.21E-05 


2.65E+001 


4.41E-01 


2.39E-03 


3.85 


2400 


q 

y 


1 


393 










3.55 


2600 


a 
o 


90 


221 


- 


- 


- 


- 


3.55 


2600 


a 
o 


SO 


221 


7.04E-07 


4.55E+000 


2.27E-01 


6.15E-04 


3.55 


2600 


a 
o 


80 


221 


1.99E-06 


1.95E+001 


3.06E-01 


1.65E-03 


3.55 


2600 


q 

y 


1 

. iU 


221 


2.96E-06 


3.68E+001 


4.69E-01 


5.06E-03 


3.70 


2600 


a 
o 


90 
.zu 


295 


1.44E-06 


1.38E+000 


3.81E-01 


5.18E-04 


3.70 


2600 


a 
o 


SO 

. JU 


295 


2.32E-06 


8.36E+000 


2.04E-01 


5.53E-04 


3.70 


2600 


a 
o 


80 
. ou 


295 


4.40E-06 


2.20E+001 


3.10E-01 


1.68E-03 


3.70 


2600 


q 

y 


1 

. iU 


295 


5.68E-06 


3.84E+001 


5.05E-01 


5.45E-03 


3.85 


2600 


a 
o 


70 
.zu 


393 


- 


- 


- 


- 


3.85 


2600 


a 
o 


SO 

. JU 


393 


6.00E-06 


1.16E+001 


2.14E-01 


5.80E-04 


3.85 


2600 


a 
o 


80 
. ou 


393 


*** 


*** 


*** 


*** 


3.85 


2600 


q 

y 


1 

. iU 


393 










3.55 


2800 


a 
o 


90 
.zu 


221 


- 


- 


- 


- 


3.55 


2800 


a 
o 


SO 

. JU 


221 


4.11E-07 


6.19E+000 


2.08E-01 


5.64E-04 


3.55 


2800 


a 
o 


80 
. ou 


221 


1.14E-06 


2.27E+001 


2.94E-01 


1.59E-03 


3.55 


2800 


q 

y 


1 

. iU 


221 


1.94E-06 


4.13E+001 


4.93E-01 


5.32E-03 


3.70 


2800 


a 
o 


70 
.zu 


295 


- 


- 


- 


- 


3.70 


2800 


a 
o 


SO 

. JU 


295 


1.91E-06 


1.13E+001 


2.21E-01 


5.99E-04 


3.70 


2800 


a 
o 


80 
. ou 


295 


2.32E-06 


2.57E+001 


2.59E-01 


1.40E-03 


3.70 


2800 


q 

y 


1 

. iU 


295 


3.78E-06 


4.24E+001 


4.99E-01 


5.38E-03 


3.85 


2800 


a 
o 


70 
.zu 


393 


- 


- 


- 


- 


3.85 


2800 


a 
o 


SO 

. JU 


393 


3.79E-06 


1.13E+001 


1.67E-01 


4.53E-04 


3.85 


2800 


a 
o 


80 
. ou 


393 


5.48E-06 


2.85E+001 


3.37E-01 


1.82E-03 


3.85 


2800 


q 

y 


1 

. iU 


393 










3.55 


3000 


a 
o 


90 
.zu 


221 


- 


- 


- 


- 


3.55 


3000 


a 
o 


SO 

. JU 


221 


4.48E-07 


8.23E+000 


2.10E-01 


5.69E-04 


3.55 


3000 


a 
o 


80 
. ou 


221 


1.05E-06 


2.49E+001 


3.66E-01 


1.98E-03 


3.55 


3000 


q 

y 


1 

. iU 


221 


2.10E-06 


3.84E+001 


5.10E-01 


5.50E-03 


3.70 


3000 


a 
o 


70 
.zu 


295 


- 


- 


- 


- 


3.70 


3000 


a 
o 


SO 

. JU 


295 


4.19E-07 


8.16E+000 


1.55E-01 


4.20E-04 


3.70 


3000 


a 
o 


80 
. ou 


295 


1.45E-06 


2.43E+001 


2.81E-01 


1.52E-03 


3.70 


3000 


q 

y 


1 

. iU 


295 


2.46E-06 


4.77E+001 


5.05E-01 


5.45E-03 


3.85 


3000 


a 
o 


90 
.zu 


393 


- 


- 


- 


- 


3.85 


3000 


a 
o 


SO 

. JU 


393 


1.80E-06 


1.31E+001 


1.56E-01 


4.23E-04 


3.85 


3000 


a 
o 


80 
. ou 


393 


3.68E-06 


3.04E+001 


3.27E-01 


1.77E-03 


3.85 


3000 


q 

y 


1 

. iU 


393 




*** 


*** 




3.55 


3200 


a 
o 


70 
.zu 


221 


- 


- 


- 


- 


3.55 


3200 


a 
o 


SO 

. JU 


221 


- 


- 


- 


- 


3.55 


3200 


a 
o 


ao 

.oU 


221 


3.63E-07 


2.77E+001 


2.58E-01 


1.40E-03 


3.55 


3200 


9 


.10 


221 


1.03E-06 


3.62E+001 


3.98E-01 


4.29E-03 


3.70 


3200 


8 


.20 


295 










3.70 


3200 


8 


.50 


295 


9.51E-07 


8.85E+000 


1.48E-01 


4.01E-04 


3.70 


3200 


8 


.80 


295 


1.64E-06 


2.56E+001 


2.01E-01 


1.09E-03 


3.70 


3200 


9 


.10 


295 


3.11E-06 


4.48E+001 


3.79E-01 


4.09E-03 


3.85 


3200 


8 


.20 


393 










3.85 


3200 


8 


.50 


393 










3.85 


3200 


8 


.80 


393 










3.85 


3200 


9 


.10 


393 
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Table 6. Input parameters (L*, T^ff, C/O, P) and the resulting mean mass loss rate, mean velocity at the outer boundary and mean 
degree of dust condensation at the outer boundary for a subset of models with = I.OMq and Aup = 2.0km s~'. The dust-to-gas 
mass ratio Pdust/Pgas is calculated from /c as described in Hofner & Dorfi (1997). Where "***" is given instead of a number, no 
meaningful wind properties could be derived due to numerical problems. 



log(L*) 




log(C-0)+12 


P 




<«out> 


(fc> 


(Pd/Pg> 




[K] 




[days] 


[Moyr-i] 


[kms-'] 






3.70 


2400 


8.20 


295 










3.70 


2400 


8.50 


295 


8.65E-07 


5.44E+00 


2.40E-01 


6.51E-04 


3.70 


2400 


8.80 


295 


2.32E-06 


2.10E+01 


2.91E-01 


1.57E-03 


3.70 


2400 


9.10 


295 


3.46E-06 


4.06E+01 


4.36E-01 


4.70E-03 


3.85 


2400 


8.20 


393 










3.85 


2400 


8.50 


393 


1.59E-06 


7.28E+00 


1.88E-01 


5.10E-04 


3.85 


2400 


8.80 


393 


4.39E-06 


2.41E+01 


2.95E-01 


1.60E-03 


3.85 


2400 


9.10 


393 


6.32E-06 


4.28E+01 


4.76E-01 


5.14E-03 


4.00 


2400 


8.20 


524 


1.09E-06 


1.48E+00 


2.26E-01 


3.07E-04 


4.00 


2400 


8.50 


524 


4.12E-06 


1.14E+01 


1.83E-01 


4.96E-04 


4.00 


2400 


8.80 


524 


9.53E-06 


2.55E+01 


3.08E-01 


1.67E-03 


4.00 


2400 


9.10 


524 


1.19E-05 


4.55E+01 


5.33E-01 


5.75E-03 


3.70 


2600 


8.20 


295 










3.70 


2600 


8.50 


295 










3.70 


2600 


8.80 


295 


1.45E-06 


2.86E+01 


3.59E-01 


1.94E-03 


3.70 


2600 


9.10 


295 


2.91E-06 


4.15E+01 


4.55E-01 


4.91E-03 


3.85 


2600 


8.20 


393 










3.85 


2600 


8.50 


393 


1.35E-06 


1.25E+01 


1.96E-01 


5.31E-04 


3.85 


2600 


8.80 


393 


2.41E-06 


2.66E+01 


2.69E-01 


1.45E-03 


3.85 


2600 


9.10 


393 


5.46E-06 


4.50E+01 


4.77E-01 


5.15E-03 


4.00 


2600 


8.20 


524 


6.72E-07 


2.27E+00 


1.99E-01 


2.70E-04 


4.00 


2600 


8.50 


524 


3.59E-06 


1.44E+01 


2.00E-01 


5.42E-04 


4.00 


2600 


8.80 


524 


5.79E-06 


2.89E+01 


2.81E-01 


1.52E-03 


4.00 


2600 


9.10 


524 


7.00E-06 


4.96E+01 


4.48E-01 


4.83E-03 


3.70 


2800 


8.20 


295 










3.70 


2800 


8.50 


295 










3.70 


2800 


8.80 


295 










3.70 


2800 


9.10 


295 


3.27E-06 


3.90E+01 


2.71E-01 


2.92E-03 


3.85 


2800 


8.20 


393 










3.85 


2800 


8.50 


393 










3.85 


2800 


8.80 


393 


9.88E-06 


3.11E+01 


2.64E-01 


1.43E-03 


3.85 


2800 


9.10 


393 


1.64E-06 


3.99E+01 


2.59E-01 


2.79E-03 


4.00 


2800 


8.20 


524 










4.00 


2800 


8.50 


524 










4.00 


2800 


8.80 


524 


2.39E-06 


2.94E+01 


2.38E-01 


1.29E-03 


4.00 


2800 


9.10 


524 


4.78E-06 


4.77E+01 


3.18E-01 


3.43E-03 


3.70 


3000 


8.20 


295 










3.70 


3000 


8.50 


295 










3.70 


3000 


8.80 


295 










3.70 


3000 


9.10 


295 










3.85 


3000 


8.20 


393 










3.85 


3000 


8.50 


393 










3.85 


3000 


8.80 


393 










3.85 


3000 


9.10 


393 










4.00 


3000 


8.20 


524 










4.00 


3000 


8.50 


524 










4.00 


3000 


8.80 


524 










4.00 


3000 


9.10 


524 


1.12E-06 


4.75E+01 


2.75E-01 


2.97E-03 


3.70 


3200 


8.20 


295 










3.70 


3200 


8.50 


295 










3.70 


3200 


8.80 


295 










3.70 


3200 


9.10 


295 










3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 










3.85 


3200 


9.10 


393 










4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 










4.00 


3200 


9.10 


524 
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Table 7. Same as|6] but for Awp = 4.0km s"' . 



log(L,) 




log(C-0)+12 P {M) <Mo„,> 


<fc> 


(Pd/Pg) 




[K] 


[days] [Moyr"'] [kms"'] 







3.70 


O y( AA 

2400 


8.20 


295 


3.70 


1 y( AA 

2400 


8.50 


295 


3.70 


2400 


8.80 


295 


3.70 


1/1 AA 

2400 


9.10 


295 


3.0J 


24U0 


O TA 
0.2U 


393 


■5 OC 

3.o5 


2400 


8.50 


393 


3.03 


'1/1 AA 

24U0 


O O A 

o.oO 


393 


o oc 

3.o5 


1y1 AA 

2400 


9.10 


393 


4.U0 


T/1 AA 


O TA 
0.2U 


J 24 


4.00 


2400 


8.50 


524 


/I AA 

4.U0 


T /I AA 


O O A 

O.oO 


524 


4.00 


2400 


9.10 


524 


3. /O 


2oU0 


O TA 

o.zU 


2^5 


3.70 


2600 


8.50 


295 


O TA 
3. /O 


2oU0 


O O A 

O.oO 


one 
2^3 


O 'TA 

3.70 


2600 


9.10 


295 


o oc 

3.o5 


2600 


8.20 


393 


■5 OC 

3.o5 


2600 


8.50 


393 


■5 OC 

3.o5 


2600 


8.80 


393 


O OC 

3.85 


2600 


9.10 


393 


A AA 

4.00 


2600 


8.20 


524 


yl AA 

4.00 


2600 


8.50 


524 


4.00 


2600 


8.80 


524 


A AA 

4.00 


2600 


9.10 


524 


1 'TA 

3.70 


10AA 

2800 


8.20 


295 


3.70 


10AA 

2800 


8.50 


295 


O 'TA 

3.70 


10AA 

2800 


8.80 


295 


1 'TA 

3.70 


10AA 

2800 


9.10 


295 


O OC 

3.85 


10AA 

2800 


8.20 


393 


■5 OC 

3.85 


10AA 

2800 


8.50 


393 


■5 OC 

3.85 


10AA 

2800 


8.80 


393 


O OC 

3.85 


10AA 

2800 


9.10 


393 


4.00 


2800 


8.20 


524 


4.00 


2800 


8.50 


524 


4.00 


2800 


8.80 


524 


A f\r\ 

4.00 


2800 


9.10 


524 


3J\J 


O AAA 

3000 


8.20 


295 


3.70 


O AAA 

3000 


8.50 


295 


■5 ^A 

3.70 


'5 AAA 

3000 


8.80 


295 


O 'TA 

3.70 


O AAA 

3000 


9.10 


295 


O OC 

3.85 


1 AAA 

3000 


8.20 


393 


■5 OC 

3.85 


'5 AAA 

3000 


8.50 


393 


■5 OC 

3.85 


'5 AAA 

3000 


8.80 


393 


O OC 

3.85 


O AAA 

3000 


9.10 


393 


A AA 

4.00 


1 AAA 

3000 


8.20 


524 


4.00 


3000 


8.50 


524 


4.00 


3000 


8.80 


524 


4.00 


3000 


9.10 


524 


3.70 


3200 


8.20 


295 


3.70 


3200 


8.50 


295 


3.70 


3200 


8.80 


295 


3.70 


3200 


9.10 


295 


3.85 


3200 


8.20 


393 


3.85 


3200 


8.50 


393 


3.85 


3200 


8.80 


393 


3.85 


3200 


9.10 


393 


4.00 


3200 


8.20 


524 


4.00 


3200 


8.50 


524 


4.00 


3200 


8.80 


524 


4.00 


3200 


9.10 


524 



1 


.05E- 


■06 


3.80E+00 


2, 


.44E- 


■01 


6 


.61E- 


■04 


3 


.02E- 


■06 


1 AAC 1 A 1 


3 


.19E- 


■01 


1 


.73E- 


■03 


4, 


.46E- 


■06 


3.86E+01 


4, 


.75E- 


■01 


5, 


.13E- 


■03 


3 


- 

.15E- 


■06 


"7 1 AA 

/. /6b+00 


2, 


- 

.13E- 


■01 


5 


- 

.77E- 


■04 


5 


.40E- 


■06 


1 1 AC 1 A 1 

2.19b+01 


2, 


.85E- 


■01 


1 


.54E- 


■03 


7 


.95E- 


■06 


A 1 CC 1 A 1 

4.15b+01 


5 


.24E- 


■01 


5 


.65E- 


■03 


2, 


.42E- 


■06 


O /;"7Tr Al 

6.6/b-Ol 


2, 


.54E- 


■01 


3 


.45E- 


■04 


8 


.50E- 


■06 


1 1 /I C 1 A 1 

1.14b+01 


2, 


.15E- 


■01 


5 


.83E- 


■04 


1 


.23E- 


■05 


1 /IOC 1 A 1 

2.42b+01 


3 


.36E- 


■01 


1 


.82E- 


■03 


1 


.59E- 


■05 


4.43E+01 


5, 


.88E- 


■01 


6 


.34E- 


■03 


7, 


.60E- 


■07 


6.24E+00 


2, 


.HE- 


■01 


5, 


.72E- 


■04 


1 


.95E- 


■06 


2.6Bb+Ul 


4, 


.OOE- 


■01 


2, 


.16E- 


■03 


3, 


.96E- 


■06 


4.39E+01 


5, 


.84E- 


■01 


6 


.30E- 


■03 


1 


- 

.67E- 


■06 


^ O /I C 1 A A 

6.24b+00 


1 


- 

.71E- 


■01 


4, 


- 

.63E- 


■04 


4, 


.04E- 


■06 


2.65E+01 


3, 


.22E- 


■01 


1 


.74E- 


■03 


5, 


.66E- 


■06 


4.65E+01 


5, 


.25E- 


■01 


5, 


.67E- 


■03 


1 


.22E- 


■06 


1.34E+00 


2, 


.06E- 


■01 


2, 


.80E- 


■04 


4, 


.85E- 


■06 


1.04E+01 


1 


.57E- 


■01 


4, 


.26E- 


■04 


8, 


.43E- 


■06 


2.80E+01 


3, 


.25E- 


■01 


1 


.76E- 


■03 


1 


.09E- 


■05 


4.55E+01 


4, 


.88E- 


■01 


5, 


.27E- 


■03 


1 


.46E- 


■06 


3.00E+01 


4, 


.55E- 


■01 


2, 


46E- 


■03 


2, 


.52E- 


■06 


4.10E+01 


4, 


.69E- 


■01 


5, 


.06E- 


■03 


1 


.49E- 


■06 


1.71E+01 


2, 


.25E- 


■01 


6, 


.lOE- 


■04 


2, 


.71E- 


■06 


2.83E+01 


3, 


.15E- 


■01 


1 


.70E- 


■03 


4, 


.77E- 


■06 


4.68E+01 


5, 


.14E- 


■01 


5, 


.55E- 


■03 


3, 


- 

.72E- 


■06 


- 

1.31E-f01 


1 


- 

.72E- 


■01 


4, 


- 

.66E- 


■04 


5 


.76E- 


■06 


2.81E+01 


2, 


.94E- 


■01 


1 


.59E- 


■03 


8, 


.07E- 


■06 


4.96E+01 


3, 


.87E- 


■01 


4, 


.18E- 


■03 


5, 


.50E- 


■07 


3.09E+01 


2, 


.84E- 


■01 


3, 


.06E- 


■03 


1 


.47E- 


■06 


3.12E+01 


3, 


.76E- 


■01 


2, 


.03E- 


■03 


3, 


.62E- 


■06 


4.27E+01 


4, 


.25E- 


■01 


4, 


.59E- 


■03 


1 


.72E- 


■06 


1.83E+01 


1 


.75E- 


■01 


4, 


.74E- 


■04 


3, 


.85E- 


■06 


2.73E+01 


2, 


.86E- 


■01 


1 


.55E- 


■03 


8, 


.OOE- 


■06 


4.76E+01 


5, 


.56E- 


■01 


6 


.OOE- 


■03 


1 


.13E- 


■06 


4.56E+01 


3, 


.22E- 


■01 


3, 


.47E- 


■03 


5, 


.17E- 


■07 


3.98E+01 


2, 


.77E- 


■01 


2, 


.99E- 


■03 


1 


.64E- 


■06 


3.33E+01 


3, 


.09E- 


■01 


1 


.67E- 


■03 


3, 


.68E- 


■06 


5.10E+01 


3, 


.65E- 


■01 


3, 


.94E- 


■03 
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Table 8. Same as|6] but for Awp = 6.0km s"' . 



log(L*) 




log(C-0)+12 


P 


{M) 


<"oiit> 


<fc> 


(Pd/Pg) 


[io] 


[K] 




[days] 


[Moyr-'] 


[km s-'] 






3.70 


2400 


8.20 


295 


- 


- 


- 


- 


3.70 


2400 


8.50 


295 


1.63E-06 


1.45E+04 


2.55E-01 


6.91E-04 


3.70 


2400 


8.80 


295 


3.81E-06 


1.92E+01 


3.22E-01 


1.74E-03 


3.70 


2400 


9.10 


295 


6.03E-06 


3.71E+01 


5.25E-01 


5.67E-03 


3.85 


2400 


8.20 


393 


- 


- 


- 


- 


3.85 


2400 


8.50 


393 


5.04E-06 


7.69E+00 


2.08E-01 


5.64E-04 


3.85 


2400 


8.80 


393 


7.81E-06 


2.68E+01 


4.47E-01 


2.42E-03 


3.85 


2400 


9.10 


393 


1.17E-05 


4.17E+01 


5.45E-01 


5.88E-03 


4.00 


2400 


8.20 


524 


- 


- 


- 


- 


4.00 


2400 


8.50 


524 


1.37E-05 


1.55E+01 


2.66E-01 


7.21E-04 


4.00 


2400 


8.80 


524 


1.60E-05 


2.76E+01 


4.94E-01 


2.67E-03 


4.00 


2400 


9.10 


524 


2.10E-05 


4.01E+01 


7.19E-01 


7.76E-03 


3.70 


2600 


8.20 


295 


- 


- 


- 


- 


3.70 


2600 


8.50 


295 


1.68E-06 


8.16E+00 


2.46E-01 


6.67E-04 


3.70 


2600 


8.80 


295 


2.93E-06 


2.42E+01 


4.03E-01 


2.18E-03 


3.70 


2600 


9.10 


295 


3.88E-06 


4.27E+01 


5.54E-01 


5.98E-03 


3.85 


2600 


8.20 


393 


- 


- 


- 


- 


3.85 


2600 


8.50 


393 


3.43E-06 


1.03E+01 


2.20E-01 


5.96E-04 


3.85 


2600 


8.80 


393 


5.48E-06 


2.79E+01 


4.16E-01 


2.25E-03 


3.85 


2600 


9.10 


393 


8.52E-06 


4.10E+01 


5.37E-01 


5.79E-03 


4.00 


2600 


8.20 


524 


2.95E-06 


1.51E+00 


2.26E-01 


3.07E-04 


4.00 


2600 


8.50 


524 


7.23E-06 


1.12E+01 


1.77E-01 


4.80E-04 


4.00 


2600 


8.80 


524 


1.07E-05 


2.84E+01 


3.67E-01 


1.98E-03 


4.00 


2600 


9.10 


524 


1.59E-05 


4.29E+01 


5.65E-01 


6.10E-03 


3.70 


2800 


8.20 


295 


- 


- 


- 


- 


3.70 


2800 


8.50 


295 


1.54E-06 


1.06E+01 


2.83E-01 


7.67E-04 


3.70 


2800 


8.80 


295 


3.12E-06 


2.50E+01 


4.33E-01 


2.34E-03 


3.70 


2800 


9.10 


295 


4.26E-06 


4.03E+01 


5.60E-01 


6.04E-03 


3.85 


2800 


8.20 


393 


- 


- 


- 


- 


3.85 


2800 


8.50 


393 


2.07E-06 


1.44E+01 


2.37E-01 


6.42E-04 


3.85 


2800 


8.80 


393 


3.74E-06 


2.71E+01 


3.15E-01 


1.70E-03 


3.85 


2800 


9.10 


393 


7.15E-06 


4.43E+01 


5.82E-01 


6.28E-03 


4.00 


2800 


8.20 


524 


- 


- 


- 


- 


4.00 


2800 


8.50 


524 


5.42E-06 


1.46E+01 


2.14E-01 


5.80E-04 


4.00 


2800 


8.80 


524 


8.76E-06 


3.10E+01 


4.11E-01 


2.22E-03 


4.00 


2800 


9.10 


524 










3.70 


3000 


8.20 


295 


- 


- 


- 


- 


3.70 


3000 


8.50 


295 


- 


- 


- 


- 


3.70 


3000 


8.80 


295 


1.16E-06 


3.13E+01 


4.49E-01 


2.43E-03 


3.70 


3000 


9.10 


295 


2.30E-06 


3.73E+01 


4.36E-01 


4.70E-03 


3.85 


3000 


8.20 


393 


- 


- 


- 


- 


3.85 


3000 


8.50 


393 


2.54E-06 


2.04E+01 


3.09E-01 


8.38E-04 


3.85 


3000 


8.80 


393 


*** 


*** 


*** 


*** 


3.85 


3000 


9.10 


393 


6.05E-06 


4.31E+01 


5.76E-01 


6.22E-03 


4.00 


3000 


8.20 


524 


- 


- 


- 


- 


4.00 


3000 


8.50 


524 


2.86E-06 


1.98E+01 


2.48E-01 


6.72E-04 


4.00 


3000 


8.80 


524 


6.37E-06 


2.71E+01 


3.21E-01 


1.74E-03 


4.00 


3000 


9.10 


524 










3.70 


3200 


8.20 


295 


- 


- 


- 


- 


3.70 


3200 


8.50 


295 


- 


- 


- 


- 


3.70 


3200 


8.80 


295 


- 


- 


- 


- 


3.70 


3200 


9.10 


295 


3.14E-06 


4.18E+01 


4.00E-01 


4.32E-03 


3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 


1.31E-06 


3.36E+01 


3.81E-01 


2.06E-03 


3.85 


3200 


9.10 


393 


3.55E-06 


4.27E+01 


4.80E-01 


5.18E-03 


4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 


2.59E-06 


9.87E+00 


3.97E-01 


1.08E-03 


4.00 


3200 


8.80 


524 


4.08E-06 


2.34E+01 


3.39E-01 


1.83E-03 


4.00 


3200 


9.10 


524 
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Table 9. Input parameters (L*, T^ff, C/O, P) and the resulting mean mass loss rate, mean velocity at the outer boundary and mean 
degree of dust condensation at the outer boundary for a subset of models with = 1.5Mq and Aup = 2.0km s~'. The dust-to-gas 
mass ratio Pdust/Pgas is calculated from /c as described in Hofner & Dorfi (1997). Where "***" is given instead of a number, no 
meaningful wind properties could be derived due to numerical problems. 



log(L*) 


^eff 


log(C-0)+12 


P 




("out) 


(fc> 


(Pd/Pg> 




[K] 




[days] 


[Moyr-i] 


[kms-'] 






3.85 


2400 


8.20 


393 










3.85 


2400 


8.50 


393 


5.40E-07 


7.95E+00 


2.53E-01 


6.86E-04 


3.85 


2400 


8.80 


393 


2.06E-06 


2.60E+01 


2.96E-01 


1.60E-03 


3.85 


2400 


9.10 


393 


4.07E-06 


4.64E+01 


4.92E-01 


5.31E-03 


4.00 


2400 


8.20 


524 










4.00 


2400 


8.50 


524 


2.23E-06 


1.17E+01 


2.30E-01 


6.23E-04 


4.00 


2400 


8.80 


524 


4.95E-06 


2.96E+01 


3.59E-01 


1.94E-03 


4.00 


2400 


9.10 


524 


8.75E-06 


4.73E+01 


5.18E-01 


5.59E-03 


4.15 


2400 


8.20 


699 










4.15 


2400 


8.50 


699 


5.26E-06 


1.83E+01 


2.68E-01 


7.26E-04 


4.15 


2400 


8.80 


699 


9.90E-06 


3.43E+01 


4.04E-01 


2.18E-03 


4.15 


2400 


9.10 


699 


1.50E-05 


5.11E+01 


5.55E-01 


5.99E-03 


3.85 


2600 


8.20 


393 










3.85 


2600 


8.50 


393 










3.85 


2600 


8.80 


393 










3.85 


2600 


9.10 


393 


*** 


*** 


*** 


*** 


4.00 


2600 


8.20 


524 










4.00 


2600 


8.50 


524 










4.00 


2600 


8.80 


524 


1.86E-06 


3.51E+01 


2.96E-01 


1.60E-03 


4.00 


2600 


9.10 


524 


2.47E-06 


5.48E+01 


3.77E-01 


4.07E-03 


4.15 


2600 


8.20 


699 










4.15 


2600 


8.50 


699 


1.50E-06 


1.63E+01 


1.89E-01 


5.12E-04 


4.15 


2600 


8.80 


699 


4.70E-06 


3.29E+01 


2.64E-01 


1.43E-03 


4.15 


2600 


9.10 


699 


5.04E-06 


5.26E+01 


3.39E-01 


3.66E-03 


3.85 


2800 


8.20 


393 










3.85 


2800 


8.50 


393 










3.85 


2800 


8.80 


393 










3.85 


2800 


9.10 


393 










4.00 


2800 


8.20 


524 










4.00 


2800 


8.50 


524 










4.00 


2800 


8.80 


524 










4.00 


2800 


9.10 


524 


9.63E-07 


4.75E+01 


2.64E-01 


2.85E-03 


4.15 


2800 


8.20 


699 










4.15 


2800 


8.50 


699 










4.15 


2800 


8.80 


699 


1.64E-06 


3.77E+01 


2.70E-01 


1.46E-03 


4.15 


2800 


9.10 


699 


2.24E-06 


5.47E+01 


2.86E-01 


3.09E-03 


3.85 


3000 


8.20 


393 










3.85 


3000 


8.50 


393 










3.85 


3000 


8.80 


393 










3.85 


3000 


9.10 


393 










4.00 


3000 


8.20 


524 










4.00 


3000 


8.50 


524 










4.00 


3000 


8.80 


524 










4.00 


3000 


9.10 


524 










4.15 


3000 


8.20 


699 










4.15 


3000 


8.50 


699 










4.15 


3000 


8.80 


699 










4.15 


3000 


9.10 


699 


7.98E-07 


4.62E+01 


1.98E-01 


2.14E-03 


3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 










3.85 


3200 


9.10 


393 










4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 










4.00 


3200 


9.10 


524 










4.15 


3200 


8.20 


699 










4.15 


3200 


8.50 


699 










4.15 


3200 


8.80 


699 










4.15 


3200 


9.10 


699 
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Table 10. Same as|9l but for Awp = 4.0km s"'. 



log(L*) 




log(C-0)+12 


P 


{M) 


<"oiit> 


<fc> 


(Pd/Pg) 


[io] 


[K] 




[days] 


[Moyr-'] 


[km s-'] 






3.85 


2400 


8.20 


393 


- 


- 


- 


- 


3.85 


2400 


8.50 


393 


1.52E-06 


6.57E+00 


2.81E-01 


7.62E-04 


3.85 


2400 


8.80 


393 


3.67E-06 


2.15E+01 


3.09E-01 


1.67E-03 


3.85 


2400 


9.10 


393 


7.23E-06 


4.64E+01 


6.87E-01 


7.41E-03 


4.00 


2400 


8.20 


524 


- 


- 


- 


- 


4.00 


2400 


8.50 


524 


3.41E-06 


6.00E+00 


2.06E-01 


5.58E-04 


4.00 


2400 


8.80 


524 


8.50E-06 


2.52E+01 


3.55E-01 


1.92E-03 


4.00 


2400 


9.10 


524 


1.22E-05 


4.76E+01 


6.56E-01 


7.08E-03 


4.15 


2400 


8.20 


699 


- 


- 


- 


- 


4.15 


2400 


8.50 


699 


8.91E-06 


l.OOE+01 


1.99E-01 


5.39E-04 


4.15 


2400 


8.80 


699 


1.72E-05 


2.93E+01 


4.04E-01 


2.18E-03 


4.15 


2400 


9.10 


699 


2.00E-05 


4.86E+01 


5.98E-01 


6.45E-03 


3.85 


2600 


8.20 


393 


- 


- 


- 


- 


3.85 


2600 


8.50 


393 


- 


- 


- 


- 


3.85 


2600 


8.80 


393 


1.71E-06 


3.51E+01 


4.16E-01 


2.25E-03 


3.85 


2600 


9.10 


393 


2.96E-06 


4.48E+01 


4.11E-01 


4.44E-03 


4.00 


2600 


8.20 


524 


- 


- 


- 


- 


4.00 


2600 


8.50 


524 


2.24E-06 


1.53E+01 


2.53E-01 


6.86E-04 


4.00 


2600 


8.80 


524 


3.87E-06 


3.13E+01 


3.78E-01 


2.04E-03 


4.00 


2600 


9.10 


524 


9.54E-06 


4.50E+01 


4.71E-01 


5.08E-03 


4.15 


2600 


8.20 


699 


- 


- 


- 


- 


4.15 


2600 


8.50 


699 


6.58E-06 


1.74E+01 


2.73E-01 


7.40E-04 


4.15 


2600 


8.80 


699 


1.21E-05 


2.79E+01 


2.95E-01 


1.60E-03 


4.15 


2600 


9.10 


699 


1.62E-05 


4.82E+01 


4.91E-01 


5.30E-03 


3.85 


2800 


8.20 


393 


- 


- 


- 


- 


3.85 


2800 


8.50 


393 


- 


- 


- 


- 


3.85 


2800 


8.80 


393 


- 


- 


- 


- 


3.85 


2800 


9.10 


393 


7.36E-07 


5.41E+01 


3.85E-01 


4.15E-03 


4.00 


2800 


8.20 


524 


- 


- 


- 


- 


4.00 


2800 


8.50 


524 


- 


- 


- 


- 


4.00 


2800 


8.80 


524 


1.77E-06 


3.87E+01 


3.56E-01 


1.93E-03 


4.00 


2800 


9.10 


524 


2.43E-06 


5.00E+01 


3.66E-01 


3.95E-03 


4.15 


2800 


8.20 


699 


- 


- 


- 


- 


4.15 


2800 


8.50 


699 


2.69E-06 


1.50E+01 


1.63E-01 


4.42E-04 


4.15 


2800 


8.80 


699 


4.29E-06 


3.58E+01 


3.87E-01 


2.09E-03 


4.15 


2800 


9.10 


699 


8.54E-06 


5.36E+01 


4.09E-01 


4.41E-03 


3.85 


3000 


8.20 


393 


- 


- 


- 


- 


3.85 


3000 


8.50 


393 


- 


- 


- 


- 


3.85 


3000 


8.80 


393 


- 


- 


- 


- 


3.85 


3000 


9.10 


393 


- 


- 


- 


- 


4.00 


3000 


8.20 


524 


- 


- 


- 


- 


4.00 


3000 


8.50 


524 


- 


- 


- 


- 


4.00 


3000 


8.80 


524 


- 


- 


- 


- 


4.00 


3000 


9.10 


524 


4.51E-08 


3.81E+01 


1.62E-01 


1.75E-03 


4.15 


3000 


8.20 


699 


- 


- 


- 


- 


4.15 


3000 


8.50 


699 


- 


- 


- 


- 


4.15 


3000 


8.80 


699 


1.02E-06 


3.70E+01 


2.34E-01 


1.27E-03 


4.15 


3000 


9.10 


699 


1.79E-06 


5.54E+01 


3.19E-01 


3.44E-03 


3.85 


3200 


8.20 


393 


- 


- 


- 


- 


3.85 


3200 


8.50 


393 


- 


- 


- 


- 


3.85 


3200 


8.80 


393 


- 


- 


- 


- 


3.85 


3200 


9.10 


393 


- 


- 


- 


- 


4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 










4.00 


3200 


9.10 


524 










4.15 


3200 


8.20 


699 










4.15 


3200 


8.50 


699 










4.15 


3200 


8.80 


699 










4.15 


3200 


9.10 


699 


9.22E-07 


4.86E+01 


2.60E-01 


2.81E-03 
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Table 11. Same as|9l but for = 6.0km s"'. 



log(L*) 




log(C-0)+12 


P 


{M) 


<"0Llt> 


<fc> 


(Pd/Ps) 


[io] 


[K] 




[days] 


[Moyr-'] 


[km s-'] 






3.85 


2400 


8.20 


393 


1.25E-06 


1.23E-01 


7.61E-01 


1.03E-03 


3.85 


2400 


8.50 


393 


4.52E-06 


1.51E+01 


3.93E-01 


1.07E-03 


3.85 


2400 


8.80 


393 


5.01E-06 


2.76E+01 


4.65E-01 


2.51E-03 


3.85 


2400 


9.10 


393 


9.82E-06 


3.92E+01 


7.73E-01 


8.34E-03 


4.00 


2400 


8.20 


524 


- 


- 


- 


- 


4.00 


2400 


8.50 


524 


7.15E-06 


1.64E+01 


3.68E-01 


9.97E-04 


4.00 


2400 


8.80 


524 


l.llE-05 


3.05E+01 


5.17E-01 


2.80E-03 


4.00 


2400 


9.10 


524 


1.47E-05 


4.55E+01 


6.10E-01 


6.58E-03 


4.15 


2400 


8.20 


699 


2.54E-05 


8.25E+00 


3.03E-01 


4.12E-04 


4.15 


2400 


8.50 


699 


2.29E-05 


1.73E+01 


4.16E-01 


1.13E-03 


4.15 


2400 


8.80 


699 


2.81E-05 


3.18E+01 


6.57E-01 


3.55E-03 


4.15 


2400 


9.10 


699 


3.34E-05 


4.53E+01 


6.37E-01 


6.87E-03 


3.85 


2600 


8.20 


393 


- 


- 


- 


- 


3.85 


2600 


8.50 


393 


2.66E-06 


l.llE+01 


3.06E-01 


8.29E-04 


3.85 


2600 


8.80 


393 


4.33E-06 


2.74E+01 


5.40E-01 


2.92E-03 


3.85 


2600 


9.10 


393 


6.22E-06 


4.52E+01 


6.70E-01 


7.23E-03 


4.00 


2600 


8.20 


524 


- 


- 


- 


- 


4.00 


2600 


8.50 


524 


6.83E-06 


2.10E+01 


4.71E-01 


1.28E-03 


4.00 


2600 


8.80 


524 


1.13E-05 


2.82E+01 


4.69E-01 


2.54E-03 


4.00 


2600 


9.10 


524 


1.36E-05 


4.44E+00 


6.30E-01 


6.80E-03 


4.15 


2600 


8.20 


699 


- 


- 


- 


- 


4.15 


2600 


8.50 


699 


1.42E-05 


2.04E+01 


3.77E-01 


1.02E-03 


4.15 


2600 


8.80 


699 


*** 




*** 


*** 


4.15 


2600 


9.10 


699 


2.03E-05 


4.66E+01 


5.50E-01 


5.93E-03 


3.85 


2800 


8.20 


393 


- 


- 


- 


- 


3.85 


2800 


8.50 


393 


- 


- 


- 


- 


3.85 


2800 


8.80 


393 


1.23E-06 


3.80E+01 


4.04E-01 


2.18E-03 


3.85 


2800 


9.10 


393 


2.07E-06 


4.47E+01 


3.84E-01 


4.14E-03 


4.00 


2800 


8.20 


524 


- 


- 


- 


- 


4.00 


2800 


8.50 


524 


2.20E-06 


1.37E+01 


2.16E-01 


5.85E-04 


4.00 


2800 


8.80 


524 


4.04E-06 


3.05E+01 


4.88E-01 


2.64E-03 


4.00 


2800 


9.10 


524 


8.55E-06 


4.24E+01 


4.18E-01 


4.51E-03 


4.15 


2800 


8.20 


699 










4.15 


2800 


8.50 


699 


*** 




*** 


*** 


4.15 


2800 


8.80 


699 


8.17E-06 


3.43E+01 


5.01E-01 


2.71E-03 


4.15 


2800 


9.10 


699 










3.85 


3000 


8.20 


393 


- 


- 


- 


- 


3.85 


3000 


8.50 


393 


- 


- 


- 


- 


3.85 


3000 


8.80 


393 


- 


- 


- 


- 


3.85 


3000 


9.10 


393 


2.04E-07 


4.54E+01 


3.30E-01 


3.56E-03 


4.00 


3000 


8.20 


524 


- 


- 


- 


- 


4.00 


3000 


8.50 


524 


- 


- 


- 


- 


4.00 


3000 


8.80 


524 


1.25E-06 


3.99E+01 


3.55E-01 


1.92E-03 


4.00 


3000 


9.10 


524 


*** 








4.15 


3000 


8.20 


699 










4.15 


3000 


8.50 


699 










4.15 


3000 


8.80 


699 










4.15 


3000 


9.10 


699 










3.85 


3200 


8.20 


393 


- 


- 


- 


- 


3.85 


3200 


8.50 


393 


- 


- 


- 


- 


3.85 


3200 


8.80 


393 


- 


- 


- 


- 


3.85 


3200 


9.10 


393 










4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 


*** 








4.00 


3200 


9.10 


524 










4.15 


3200 


8.20 


699 










4.15 


3200 


8.50 


699 










4.15 


3200 


8.80 


699 










4.15 


3200 


9.10 


699 
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Table 12. Input parameters (L*, Tetf , C/O, P) and the resulting mean mass loss rate, mean velocity at the outer boundary and mean 
degree of dust condensation at the outer boundary for a subset of models with = I.OMq and Aup = 2.0km s~'. The dust-to-gas 
mass ratio Pdust/Pgas is calculated from /c as described in Hofner & Dorfi (1997). 



log(L,) 




log(C-0)+12 


P 


<M> 


<Woul> 


<fc) 


(Pd/Pg) 


[£«] 


[K] 




[days] 


[Meyr-'l 


[km s-'] 






3.85 


2400 


8.20 


393 










3.85 


2400 


8.50 


393 










3.85 


2400 


8.80 


393 










3.85 


2400 


9.10 


393 


1.57E-06 


5.38E+01 


4.34E-01 


4.68E-03 


4.00 


2400 


8.20 


524 










4.00 


2400 


8.50 


524 










4.00 


2400 


8.80 


524 


2.05E-06 


3.29E+01 


3.23E-01 


1.75E-03 


4.00 


2400 


9.10 


524 


*** 


*** 


*** 


*** 


4.15 


2400 


8.20 


699 










4.15 


2400 


8.50 


699 


*** 


*** 


*** 


*** 


4.15 


2400 


8.80 


699 


4.13E-06 


3.68E+01 


3.52E-01 


1.90E-03 


4.15 


2400 


9.10 


699 


6.13E-06 


5.59E+01 


4.28E-01 


4.62E-03 


3.85 


2600 


8.20 


393 










3.85 


2600 


8.50 


393 










3.85 


2600 


8.80 


393 










3.85 


2600 


9.10 


393 










4.00 


2600 


8.20 


524 










4.00 


2600 


8.50 


524 










4.00 


2600 


8.80 


524 










4.00 


2600 


9.10 


524 


1.27E-06 


5.09E+01 


3.08E-01 


3.32E-03 


4.15 


2600 


8.20 


699 










4.15 


2600 


8.50 


699 










4.15 


2600 


8.80 


699 










4.15 


2600 


9.10 


699 


2.33E-06 


5.29E+01 


2.83E-01 


3.05E-03 


3.85 


2800 


8.20 


393 










3.85 


2800 


8.50 


393 










3.85 


2800 


8.80 


393 










3.85 


2800 


9.10 


393 










4.00 


2800 


8.20 


524 










4.00 


2800 


8.50 


524 










4.00 


2800 


8.80 


524 










4.00 


2800 


9.10 


524 










4.15 


2800 


8.20 


699 










4.15 


2800 


8.50 


699 










4.15 


2800 


8.80 


699 










4.15 


2800 


9.10 


699 










3.85 


3000 


8.20 


393 










3.85 


3000 


8.50 


393 










3.85 


3000 


8.80 


393 










3.85 


3000 


9.10 


393 










4.00 


3000 


8.20 


524 










4.00 


3000 


8.50 


524 










4.00 


3000 


8.80 


524 










4.00 


3000 


9.10 


524 










4.15 


3000 


8.20 


699 










4.15 


3000 


8.50 


699 










4.15 


3000 


8.80 


699 










4.15 


3000 


9.10 


699 










3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 










3.85 


3200 


9.10 


393 










4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 










4.00 


3200 


9.10 


524 










4.15 


3200 


8.20 


699 










4.15 


3200 


8.50 


699 










4.15 


3200 


8.80 


699 










4.15 


3200 


9.10 


699 
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Table 13. Same as[T2l but for Awp = 4.0km s"'. 



log(L,) 




log(C-0)+12 P {M) <Mo„,> 


<fc> 


(Pd/Pg) 




[K] 


[days] [Moyr"'] [kms"'] 







i.o5 


O /I AA 

2400 


8.20 


393 










■5 OC 

i.o5 


1 /I AA 

2400 


8.50 


393 










3.85 


1 /I AA 

2400 


8.80 


393 


2.09b-06 


2.8 /b+Ul 


3.6sb-01 


1.99b-03 


O OC 

3.85 


l/l AA 

2400 


9.10 


393 


2.4sb-06 


e 1 /I c 1 A 1 

5.14b+Ul 


5.36E-01 


3. /ob-U3 


4.UU 


T/1 AA 


O TA 
0.2U 


eT A 
524 


1 "yec A/i 
1. /5rl-06 


1 AOC 1 Al 

l.Uob+Ul 


T OiiC Al 

2.o6b-Ul 


1 one c\A 


4.00 


2400 


8.50 


524 


1.68E-06 


1.13E+01 


r\ rvic Al 

2.93b-01 


/.y4b-U4 


4.U0 


T /I AA 


O O A 

o.oO 


524 


/I T 1 C A/i 

4. / lil-U6 


T TOC 1 A 1 

3.3Sb+Ul 


A A iiC A 1 

4.46b-Ul 


o /1 1 c m 
Z.41rl-Uj 


4.00 


2400 


9.10 


524 


C O /I C 

6.84b-06 


e 1 AC 1 A 1 

5.10b+01 


5.46E-01 


r one m 


/lie 

4.15 


0/1 AA 

/4U0 


O TA 


d99 


"7 AOC A/i 

/.L)21l-U6 


1 A 1 C 1 A 1 


Q TQC Al 

3. /3b-Ul 


j.U /b-U4 


yl 1 C 

4.15 


2400 


8.50 


699 


5.34E-06 


1.65E+01 


1 AAC A 1 

3.00b-01 


O 1 Ayl 


/lie 
4.15 


T /I AA 

z4U0 


O O A 

O.oO 




A TOC A/i 

y.z/ii-u6 


T A AC 1 A 1 


A AiiC A 1 

4.y6b-ui 


z.oob-U3 


4.15 


2400 


9.10 


699 


1.68E-05 


A o ec 1 A 1 

4.85b+0i 


5.85E-01 


6.31b-03 


1 OC 

3.85 


2oU0 


O TA 


393 










■5 OC 

3.85 


2600 


8.50 


393 










1 OC 

3.85 


2oU0 


O O A 

O.oO 


393 










IOC 

3.85 


2600 


9.10 


393 


rv one Csn 

9.89b-07 


e /I /I c 1 A 1 
5.44b+Ul 


A 1 AC A 1 

4.30b-01 


A ^/1C AT 

4.64b-03 


4.00 


2600 


8.20 


524 










4.00 


2600 


8.50 


524 










4.00 


2600 


8.80 


524 


1 OOC A/C 

1.82b-06 


A AOC 1 Al 

4.02b+Ul 


A AAC A 1 

4.00b-01 


T 1 /^C Al 

2.16b-03 


4.00 


2600 


9.10 


524 


2.44b-06 


e T AC 1 A 1 

5. /Ob+01 


A A t^T2 A 1 

4.46b-01 


/I O 1 C AT 

4.81b-03 


4.15 


2600 


8.20 


699 










/lie 

4.15 


2600 


8.50 


699 










4.15 


2600 


8.80 


699 


4.19E-06 


3.88E+01 


4.11E-01 


2.22E-03 


/lie 

4.15 


2600 


9.10 


699 


9.12b-06 








O OC 

3.85 


'^OAA 

2800 


8.20 


393 










■5 OC 

3.85 


/^OAA 

2800 


8.50 


393 










'5 OC 

3.85 


10AA 

2800 


8.80 


393 










IOC 

3.85 


'^OAA 

2800 


9.10 


393 










4.00 


TOAA 

2800 


8.20 


524 










4.00 


2800 


8.50 


524 










4.00 


2800 


8.80 


524 










4.00 


2800 


9.10 


524 










/lie 

4.15 


2800 


8.20 


699 










/lie 

4.15 


2800 


8.50 


699 










4.15 


2800 


8.80 


699 










/lie 

4.15 


2800 


9.10 


699 


1 TAC 

1.39b-06 


5.90b+01 


1 O AC A 1 

3.20b-01 


3.45b-03 


IOC 

3.85 


O AAA 

3000 


8.20 


393 










O OC 

3.85 


O AAA 

3000 


8.50 


393 










■5 OC 

3.85 


'5 AAA 

3000 


8.80 


393 










O OC 

3.85 


O AAA 

3000 


9.10 


393 










/I AA 

4.00 


1 AAA 

3000 


8.20 


524 










4.00 


3000 


8.50 


524 










4.00 


3000 


8.80 


524 










4.00 


3000 


9.10 


524 










4.15 


3000 


8.20 


699 










/lie 

4.15 


3000 


8.50 


699 










/lie 

4.15 


3000 


8.80 


699 










/lie 

4.15 


3000 


9.10 


699 










3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 










3.85 


3200 


9.10 


393 










4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 










4.00 


3200 


9.10 


524 










4.15 


3200 


8.20 


699 










4.15 


3200 


8.50 


699 










4.15 


3200 


8.80 


699 










4.15 


3200 


9.10 


699 
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Table 14. Same as[T2l but for Awp = 6.0km s"'. 



log(L,) 




log(C-0)+12 P {M) <Mo„,> 


<fc> 


(Pd/Ps) 




[K] 


[days] [Moyr"'] [kms"'] 







i.o5 


1 /I AA 

2400 


8.20 


393 


1.03b-06 


0. / /b+00 


3.16E-01 


A TAC r\A 


■5 OC 

i.o5 


/^/l AA 

2400 


8.50 


393 


(\ ICC A"7 

9.36h-0/ 


O IOC 1 AA 

s.28b+00 


T A/IC Al 

2.96b-01 


O AOC A/1 


3.85 


"^/l AA 

2400 


8.80 


393 


3.99b-06 


1 AAC 1 A 1 

3.U0b+Ul 


e O AC A 1 

5.oOb-01 


1 i ATH Al 

3.14b-U3 


o oc 

3.85 


l/l AA 

2400 


9.10 


393 


5.9sb-06 


A A 1 C 1 A 1 

4.Ulb+Ul 


6.60E-01 


n 1 Tc Al 

/.lzb-U3 


A AA 

4.00 


O/) AA 

2400 


8.20 


524 


/.39b-06 


1.65E+01 


T "7/lC Al 

3. /6b-01 


5.1 lE-04 


4.00 


2400 


8.50 


524 


o.3zb-Uo 


1.55E+01 


A AAC A 1 

4.00b-01 


1 AOC Al 


4.00 


2400 


8.80 


524 


9.ozb-Uo 


n AAC 1 A 1 

2.99b+01 


5.16E-01 


/-t TAC Al 

2.7yb-03 


4.00 


2400 


9.10 


524 


1 TIC AC 

1.33b-05 


A 1 ec 1 A 1 

4.15b+01 


n n cm A 1 

/. /5b-01 


O l/^C Al 

8.36b-03 


/lie 

4.15 


2400 


8.20 


699 


1 A/1 

l.z/b-06 


o Tec 1 A 1 

2.3jb+01 


5.14E-01 


/I AOC A/1 

6.95b-04 


yl 1 C 

4.15 


2400 


8.50 


699 


1.43E-05 


1 ^ /I C 1 A 1 

l.o4b+01 


5.74E-01 


1.56E-03 


4.15 


2400 


8.80 


699 


l.lOE-05 


1 A 1 C 1 A 1 

3.91b+Ul 


5.65E-01 


1 A/^C Al 

3.06b-03 


4.15 


2400 


9.10 


699 


n CQC AC 

2.53b-05 


A ^ OC I A 1 

4.Dsb+0i 


7.74b-01 


O 1 CC Al 

8.35b-03 


o oc 

3.85 


T^AA 

2600 


8.20 


393 










■5 oc 

3.85 


2600 


8.50 


393 










o oc 

3.85 


2600 


8.80 


393 


1.17E-06 


3.65E+01 


4.25E-01 


T 1 AC Al 

2.30b-03 


3.85 


2600 


9.10 


393 


1.58E-06 


5.54E+01 


5.10E-01 


5.50E-03 


4.00 


2600 


8.20 


524 










4.00 


2600 


8.50 


524 










4.00 


2600 


8.80 


524 


3.08b-06 


1 ACC r A 1 

3.95b+01 


A A'TC A 1 

4.9/b-Ol 


/-t zIAC Al 

2.69b-03 


4.00 


2600 


9.10 


524 


A A/1 


A O 1 C r A 1 


4.56E-01 


A ATC Al 

4.yzb-U3 


yl 1 C 

4.15 


2600 


8.20 


699 


7.65E-06 


i-<i T OC 1 A 1 


A T /I C A 1 

4.34b-01 


C AAC f\A 


4.15 


2600 


8.50 


699 


6.65E-06 


2.12E+01 


3.61E-01 


9.78E-04 


/lie 

4.15 


2600 


8.80 


699 


8.09b-06 


1 C AC 1 A 1 

3.50b+01 


4.26E-01 


/-t 1 AC Al 

2.30b-03 


/lie 

4.15 


2600 


9.10 


699 


1 .42E-05 


A AOC 1 A 1 

4.9/b+Ol 


A Aec Al 

4.9jb-01 


c n AU Al 

5.34b-03 


o oc 

3.85 


2800 


8.20 


393 










■5 OC 

3.85 


/^OAA 

2800 


8.50 


393 










o oc 

3.85 


2800 


8.80 


393 










o oc 

3.85 


'^OAA 

2800 


9.10 


393 










/I AA 

4.00 


2800 


8.20 


524 










4.00 


2800 


8.50 


524 










4.00 


2800 


8.80 


524 


T T7C A/1 

3.2 /b-Oo 


1 O AC 1 A 1 

2.o0b+01 


T Aec Al 

2.05b-01 


l.llE-03 


4.00 


2800 


9.10 


524 


1 nnu A^ 

1.77b-06 


5.57E+01 


1 O oc A 1 

3.88b-01 


A 1 AC AT 

4.iyb-03 


/lie 

4.15 


2800 


8.20 


699 










4.15 


2800 


8.50 


699 










4.15 


2800 


8.80 


699 


3.8sb-06 


A /I ec 1 A 1 

4.4Db+Ul 


A AQC Al 

4.03b-01 


1 1 oc AT 

2.l6b-03 


4.15 


2800 


9.10 


699 


c 1 AC r\c 

5.30b-06 


e mc 1 A 1 


A AAC A 1 

4.09b-01 


A /lie AT 

4.41b-03 


IOC 

3.85 


O AAA 

3000 


8.20 


393 










o oc 

3.85 


O AAA 

3000 


8.50 


393 










o oc 

3.85 


O AAA 

3000 


8.80 


393 










o oc 

3.85 


O AAA 

3000 


9.10 


393 










,1 AA 

4.00 


O AAA 

3000 


8.20 


524 










4.00 


3000 


8.50 


524 










4.00 


3000 


8.80 


524 










4.00 


3000 


9.10 


524 










4.15 


3000 


8.20 


699 










/lie 

4.15 


3000 


8.50 


699 










4.15 


3000 


8.80 


699 










4.15 


3000 


9.10 


699 


1 .76E-06 


A O AC 1 A 1 

4.s0b+0i 


/-t CTC Al 

z.o /b-Ul 


1 1 AC Al 

3.1Ub-U3 


3.85 


3200 


8.20 


393 










3.85 


3200 


8.50 


393 










3.85 


3200 


8.80 


393 










3.85 


3200 


9.10 


393 










4.00 


3200 


8.20 


524 










4.00 


3200 


8.50 


524 










4.00 


3200 


8.80 


524 










4.00 


3200 


9.10 


524 










4.15 


3200 


8.20 


699 










4.15 


3200 


8.50 


699 










4.15 


3200 


8.80 


699 










4.15 


3200 


9.10 


699 











